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GENERAL  INTRODUCTION. 


Study  of  turbulence  is  of  vital  scientific,  military  and  economic  interest.  Advances  in  several 
areas  of  external  aerodynamics  and  internal  combustor  flows  of  interest  to  Air  Force  hinge  on 
our  ability  to  clearly  understand  and  adequately  predict  the  effects  of  turbulence.  At  the  current 
time,  however,  there  exists  a  substantial  gap  between  our  knowledge  of  the  physics  of  turbulence 
phenomenon  and  the  physics  that  is  incorporated  into  turbulence  models,  especially  subgrid 
closures.  The  disciplines  of  turbulence  theoiy/analysis  (e.g.,  rapid  distortion  theory,  spectral 
closure  models),  high-order  turbulence  modeling  (e.g.,  second-moment  closures,  structure-based 
models  and  realizability  constraints)  and  turbulence  simulation  (DNS-  direct  numerical 
simulations,  and  LES-  large  eddy  simulations)  are  evolving  independently  with  very  little  cross 
fertilization  of  ideas.  For  example,  the  currently  popular  LES  subgrid  closures  (e.g., 
Smagorinsky,  dynamic  Smagorinsky)  are  algebraic  in  nature;  completely  insensitive  to  extra 
rates  of  strain  such  as  rotation,  curvature,  and  buoyancy  and,  further,  may  not  even  be  realizable. 
These  major  deficiencies  in  the  LES-SGS  modeling  are  tolerated  despite  the  fact  that,  in  higher- 
order  closures,  fiiese  physical  effects  and  mathematical  constraints  have  long  been  represented 
adequately.  Further,  we  would  like  to  point  out  that  the  very  premise  of  detached-eddy 
simulation  (DES)  approach  -  that  is  seen  as  the  practical  computational  tool  for  turbulence  -  is 
erroneous.  This  is  due  to  the  fact  that  that  inhomogeneous  spatial  filtering  i  s  inevitable  in  this 
method,  and  yet  the  governing  equations  ignore  the  effects  that  necessarily  arise  with 
inhomogeneous  filtering  of  the  velocity  field. 

Much  of  the  disparity  between  what  is  known  of  turbulence  physics  and  what  is  incorporated 
into  the  LES  models  is  due  to  the  reason  that  the  frameworks  of  various  developments  have 
traditionally  been  vastly  different.  The  objective  of  the  research  performed  under  this  grant  was 
to  find  means  of  unifying  these  disciplines  on  a  common  framework  such  that  turbulence  closure 
models  -  especially,  LES-SGS  closures  -  can  benefit  from  the  decades  of  wisdom  and  experience 
incumbent  in  turbulence  theories  and  high-order  closures. 

The  flowchart  given  below  summarizes  the  plaimed  sequence  of  investigations  that  is  purported 
to  result  in  a  turbulence  closure  model  that  can  be  used  at  any  degree  of  evolution.  We  start  with 
the  premise  that  turbulence  modeling  can  be  accomplished  with  some  degree  of  fidelity  to  file 
governing  Navier- Stokes  equations  only  at  the  second -order  closure  level.  At  this  level,  different 
turbulence  theories  and  analyses  like  RDT  (rapid  distortion  theory),  realizability  and  effects  of 
extra  rates  of  strain  (via  the  Poisson  equation  for  fluctuating  pressure)  can  formally  contribute  to 
the  closure  models  of  rapid  pressure-strain  correlation  and  dissipation.  In  the  flowchart,  the 
rectangular  boxes  in  light  blue  represent  the  fundamental  theory/analysis  studies  that  directly 
contribute  to  the  seven-equation  second-order  (RANS)  turbulence  modeling  closure  (dark-blue 
box).  Once  the  RANS  closure  has  been  achieved,  one  can  (proposed  future  work)  systematically 
develop  a  seven-equation  PANS  (partially-averaged  Navier-Stokes)  model  to  represent  the 
effects  of  partially-averaged  fields.  Under  work  funded  mostly  by  NASA,  we  have  already 
established  the  PANS  protocol  for  simple  two-equation  models.  This  experience  will  definitely 
aid  us  in  the  development  of  the  PANS  protocol  for  seven-equation  closure.  As  the  seveit 
equation  closure  for  PANS  or  LES-SGS  is  certain  to  be  infeasible,  we  will  have  to  effect 
reduction  to  a  smaller  equation  set.  This  can  be  accomplished  by  employing  algebraic  stress 
modeling  approach  fiiat  has  been  perfected  in  our  group.  This  will  ultimately  result  in  the 


derivation  of  an  algebraic  PANS  closure  model  that  can  be  used  in  the  context  of  PANS  or  even 
LES.  Thus,  the  body  of  turbulence  knowledge  developed  over  the  years  in  turbulence  theoiy- 
analysis-modeling  studies  can  be  brought  to  bear  on  computing  filtered  velocity  fields.  With  such 
models,  any  new  physical  effect  (rotation,  curvature,  buoyancy)  can  be  systematically  introduced 
into  model  closure  with  resorting  to  any  kind  of  ad  hoc  modifications. 


Unifying  Framework 


PANS 


LES  DNS 


In  the  completed  research,  three  important  topics  that  make  crucial  contributions  to  the  above 
model  development  strategy  were  investigated:  (i)  new  realizability  constraints  for  rapid 
pressure-strain  correlation  modeling;  (ii)  new  physbal  modeling  constraints  that  are  needed  for 
ensuring  consistency  with  rapid  distortion  theory;  and  (iii)  turbulent  transport  models  for 
partially-averaged  (sub-grid)  kinetic  energy  and  dissipation.  We  will  first  provide  brief 


descriptions  of  each  of  the  topic  and  present  the  detailed  reports  from  each  study  in  the  following 
chapters. 


New  realizability  constraints.  In  second- moment-closure  method  in  turbulence  modeling, 
measures  to  insure  a  realizable  turbulence  model  are  currently  limited  to  constraining  the 
Reynolds  stress  to  physically  possible  values.  These  constraints  address  neither  the  realizability 
of  the  other  modified  statistical  moments  (e.g.,  pressure-strain  correlation)  nor  the  underlying 
causes  of  unrealizable  Reynolds  stresses.  For  achieving  increased  consistency  with  flow  physics 
in  SMC,  we  propose  the  additional  requirement  that  the  closure  model  for  each  of  the  xmclosed 
statistical  moments  in  the  Reynolds  stress  equation  be  individually  realizable.  We  then  proceed 
to  derive  two  realizability  constraints  on  the  rapid  pressure-strain  term:  (i)  the  so-called  Mijkl 
tensor  must  be  positive  semi-definite  and  (ii)  die  correlation  closure  must  satisfy  Schwarz 
inequality.  Calculations  with  currently-popular  models  show  that  unrealizable  rapid  pressure- 
strain  correlation  precedes  unrealizable  Reynolds  stress.  It  is  also  demonstrated  that  when  the 
Launder,  Reece  and  Rodi  (LRR)  rapid  pressure-strain  correlation  model  is  modified  (truncated) 
to  satisfy  the  new  constraints,  Reynolds  stress  realizability  is  always  preserved.  These  findings 
clearly  indicate  that  unrealizable  closure  model  is  the  cause  of  Reynolds  stress  realization 
violation  and  highlight  the  importance  of  the  new  constraints.  Incidentally,  the  truncated  LRR 
model  -  which  is  folly  realizable  and  piece- wise  linear  -  comes  close  to  satisfying  two  crucial 
requirements  of  a  rapid  pressure-strain  closure. 


New  physical  constraints  based  on  RDT.  Rapid  distortion  calculations  of  initially  anisotropic 
turbulence  are  performed  to  better  understand  the  physics  of  the  pressure-strain  correlation  in 
strain  dominated  flows.  Based  on  the  results  of  simulations  we  infer  important  physical 

characteristics  of  the  “rapid”  pressure- strain  correlation  in  mean  strain-dominated 

turbulence:  i)  it  vanishes  when  foere  is  no  production  of  anisotropy,  ii)  it  tends  to  decrease 
Reynolds  stress  anisotropy,  and  iii)  its  magnitude  is  mostly  smaller  than  that  of  production.  The 
observed  characteristics  are  proposed  as  criteria  that  pressure-strain  correlation  models  should 

satisfy.  Although  these  properties  of  appear  intuitively  self-evident,  all  of  the  current 

popular  models  violate  foe  above  criteria  for  a  sizeable  subset  of  anisotropic  initial  conditions. 
Reynolds  stress  transport  model  calculations  show  that  unphysical  (inconsistent)  and 
unrealizable  model  behavior  can  be  directly  attributed  to  these  violations. 

Transport  models  for  subgrid  kinetic  energy  and  dissipation.  If  foe  filtered  stresses  are  to  be 
calculated  from  transport  equations  (rather  than  algebraic  closures),  processes  such  as  pressure- 
strain  correlation,  turbulent  transport  and  dissipation  (in  foe  filtered  equation  context)  need  to  be 
modeled.  In  previous ’work  (performed  under  NASA  grant),  we  have  demonstrated  that  closure 
for  all  processes,  except,  transport  can  be  inferred  from  foe  parent  RANS  model  by  simply 
invoking  foe  averaging- invariance  principal  in  a  consistent  maimer.  The  turbulent  transport 
model,  however,  requires  independent  calibration. 


Chapter  1  New  perspectives  on  realizability  of  pressure -strain 
correlation  closures 


1.1  Introduction 

Realizability  requirement  enunciates  the  rudimentary  expectation  that  an  acceptable  turbulence 
model  closure  expression  be  based  on  the  statistics  of  the  velocity  field  that  is  physically 
achievable  or  realizable.  The  covariance  tensor,  or  the  Reynolds  stress,  of  any  velocity  field 
governed  by  the  Navier-Stokes  equations  is  positive  semi-definite  exhibiting  two  important 
characters;  the  diagonal  components  (energies)  are  non-negative  and  the  off-diagonal 
components  satisfy  the  Schwarz  inequality.  For  second  moment-closures  (SMC),  the  realizability 
constraint  as  proposed  by  Schumann  (1977)  requires  that  a  turbulence  model  yield  Reynolds 
stresses  that  satisfy  these  conditions  (see  also  Lumley,  1978).  Over  the  last  two  decades, 
Schumaim  realizability  constraint  has  served  as  the  theoretical  basis  for  several  turbulence 
models  (e.g.,  Ristorcelli,  Lumley  and  Abid,  1995;  Johansson  and  Hallback,  1994;  and  Sjogren 
and  Johansson,  2000).  In  this  paper,  we  revisit  the  realizability  issue  and  present  a  new 
perspective. 

We  suggest  that  a  SMC  turbulence  model  must  not  only  yield  a  positive  semi-definite  Reynolds 
stress  tensor,  but  the  model  expressions  for  individual  unclosed  moments  (e.g.,  pressure-strain 
correlations)  must  themselves  be  realizable.  While  the  theoretical  justification  for  this  extended 
or  comprehensive  realizability  condition  is  to  achieve  improved  consistency  with  turbulence 
physics,  there  is  also  an  important  practical  consideration  In  SMC  (second-moment  closure) 
methodology,  Reynolds  stress  is  calculated  from  an  evolution  equation.  The  evolution  rate  is 
composed  of  contributions  from  production,  pressure-strain  correlation,  turbulent  transport, 
dissipation  and  other  processes.  Most  of  these  terms  require  closure  modeling.  It  is  logical  then 
fliat  the  most  probable  of  umealizable  calculated  stress  is  unphysical  models  fort  the  individual 
terms  contributing  towards  the  evolution  rate.  One  of  the  most  fundamental  prerequisites  of  a 
physically  acceptable  closure  model  is  that  it  be  realizable.  Each  unclosed  SMC  term  is  a 
statistical  moment,  the  closure  model  for  which  must  satisfy  certain  realizability  constraints.  For 
example,  an  uiurealizable  pressure-strain  correlation  model  is  one  in  which  the  closure  expression 
cannot  be  possibly  obtained  from  any  physically  permissible  (real)  pressure  and  velocity  fields. 
Closure  models  not  satisfying  the  appropriate  realizability  constraints  (given  in  Section  1 .2)  must 
be  considered  unphysical  and  unacceptable. 

The  current  approach  attempts  to  enforce  reliazability  of  Reynolds  stress  without  due 
consideration  of  the  reliazability  of  the  underlying  closure  models.  Although  such  a  model  may 
lead  to  apparently  plausible  Reynolds  stresses,  the  dynamics  of  SMC  models  will  certainly  be 
inconsistent  with  turbulence  physics.  Such  realizabilty  enforcement  will  defeat  the  ultimate  goal 
which  is  to  develop  turbulence  models  of  high  degree  of  fidelity  to  the  Navier-Stokes  equation. 

The  flow  physics  underlying  the  new  constraints  can  be  explained  as  follows.  All  velocity  fields 
that  are  solutions  of  the  Navier-Stokes  equations  are  real-valued  and,  hence,  yield  a  covariance 
matrix  (Re3molds  stress)  whose  diagonal  components  are  positive  and  off-diagonal  components 
satisfy  die  Schwarz  inequality.  However,  not  all  real-valued  fields  with  covariance  tensor 
satisfying  the  above  conditions  can  be  solutions  of  the  flow  equations.  A  permissible  flow- filed 


solution  must  be  real- valued  and  satisfy  additional  kinematic  criteria  such  as  the  appropriate 
compressibility  condition.  For  example,  in  incompressible  flow,  the  real- valued  velocity  field 
must  also  be  divergence- free.  The  new  constraints  proposed  in  this  paper  go  a  step  beyond  and 
attempt  to  eliminate  real- valued  velocity  fields  diat  do  not  satisfy  the  additional  requirements. 

In  this  paper  we  will  focus  on  the  reliazabilily  issues  of  rapid  pressure-strain  correlation  term. 
Exploiting  the  divergence- free  character  of  an  incompressible  flow  field,  we  derive  reliazability 
constraints  on  the  rapid  pressure- strain  correlation  term.  Simple  analysis  shows  that  many 
popular  linear  and  non-linear  models  do  not  satisfy  there  constraints.  Detailed  model  calculations 
are  performed  in  homogeneous  turbulence  to  demonstrate  unequivocally  that  the  onset  of 
unrealizable  Reynolds  stress  is  preceded  by  unrealizable  rapid  pressure-strain  correlation.  This 
establishes  the  cause-effect  relationship  between  umealizable  pressure-strain  correlation  model 
and  unrealizable  Reynolds  stress  highlighting  the  importance  of  the  new  constraints. 

An  introduction  to  the  rapid  pressure-strain  correlation  modeling  is  given  in  section  2.  In  section 
1.3,  we  discuss  the  current  reliazability  approach  and  motivate  die  need  for  additional 
constraints.  The  new  realizability  constraints  on  the  rapid-pressure  statistics  are  derived  in 
section  1.4.  Section  1.5  contains  the  solution  procedure  and  the  results  are  presented  in  section 
1.6.  We  close  with  a  summary  and  conclusions  section  1.7. 

1.2  Rapid  pressure-strain  correlation 
The  Reynolds  stress  evolution  is  equation  by: 


DUfUj 

Dt 


(1.2.1) 


where  u^Uj  is  the  Reynolds  stress.  The  various  terms  in  the  above  equation  are  the  time-rate-of- 
change  following  the  mean  flow,  production,  pressure- strain  correlation,  dissipation  and 
transport.  Several  of  these  require  closure  modeling. 

Pressure-strain  correlation  modeling  mediodology  for  incompressible  flows  is  derived  from  the 
early  work  of  Chou  (1945).  We  start  with  the  Poisson  equation  for  the  fluctuating  pressure: 


£p_ 

dXfX, 


d\p'  +  P^) 

dx^Xi 


=  -2(S,jSji  +  WyWj,)  -  (SySj,  +  WyWji) 


(1.2.2) 


Here,  and  throughout  the  remainder  of  the  paper,  the  following  notation  will  be  used:  Sy  and  Wy 
represent  the  strain  and  rotation  rates  of  mean  flow;  s  y  and  wy  represent  those  of  the  fluctuating 
filed.  Repeated  English  indices  will  imply  summation.  The  pressure  fluctuations  are  composed  of 
two  parts  p=p^+  p' :  the  slow  part  (p®)  is  due  to  the  fluctuating  velocity  field  alone,  and  the  rapid 
part  (p’)  is  due  to  the  interaction  between  die  mean  and  fluctuation  velocity  fields.  Accordingly, 
it  is  normal  practice  to  decompose  the  pressure- strain  correlation  into  slow  and  rapid  parts 

(0„  =  0*  +  (p’j )  and  model  each  component  separately  (Launder,  Reece  and  Rodi,  1975). 


The  Poisson  equation  for  fluctuating  rapid  pressure  can  be  written  as 


(1.2.3) 


|^=-2(v„ 

dx^x^ 

Using  the  Greens  function  solution  to  the  Poisson  equation,  it  can  be  shown  tiiat  4ie  (rapici) 
pressure  velocity-gradient  correlation  tensor  can  be  written  in  the  form 


P 


dxj 
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The  fourth-order  tensor  Miju  is  given  by 


1  fl 
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(1.2.4) 


(1.2.5) 


where  Rn  is  the  two-point  velocity  correlation:  )m,5^  -t  ) .  It  is  ttien  straight-forward 

to  show  that  the  rapid  pressure-strain  correlation  can  be  expressed  as 

=  (1.2.6) 

Xi 

Modeling  rapid  pressure-strain  correlation  is  tantamount  to  modeling  the  Mija  tensor.  Complete 
description  of  the  Mya  tensor  requires  the  knowledge  of  the  Reynolds  stress  (compatibility 
tensor)  and  the  so-called  dimensionality  tensor  (Dy)  as  is  well  explained  in  Kassinos,  Reynolds 
and  Rogers  (2001).  The  evolution  of  dimensionality  tensor  is  governed,  to  a  large  extent,  by  the 
initial  dimensionality  and  the  mean- velocity  gradients.  The  dimensionality  tensor  is  unclosed  in 
terms  of  the  single-point  tensors  used  in  the  conventional  SMC  methods.  Fort  the  given 
Reynolds  stress,  the  rapid  pressure-strain  correlation  can  take  a  range  of  values  depending  on  the 
dimensionality  of  turbulence.  This  is  the  crux  of  the  closure  problem.  While  carrying  evolution 
equations  for  dimensionality  is  an  option  it  is  not  generally  undertaken  due  to  the  added 
computational  expense.  The  practice  in  traditional  SMC  modeling  is  to  postulate  a  closure 
expression  for  the  Myu  tensor  that  is  based  on  some  implied  average  value  of  dimensionality 
tensor.  In  effect,  the  conditional  (based  on  a  specific  value  of  dimensionality  tensor)  correlation 
between  pressure  and  strain  fluctuations  is  modeled  with  the  unconditional  correlation  (based  on 
some  average  dimensionality  tensor).  Thus,  the  need  fer  dimensionality  tensor  is  obviated,  at  the 
expense  of  accuracy. 

Lacking  knowledge  of  the  dimensionality  tensor,  the  closure  expression  for  pressure- strain 
correlation  is  generally  required  to  satisfy  several  physical  and  mathematical  constraints  to 
ensure  some  degree  of  fidelity  of  the  flow  physics.  A  complete  list  of  currently  mandated 
constraints  can  be  found  in  Pope  (2000).  Three  of  the  important  constraints  are  that  rapid 
pressure-strain  correlation  model  must  (i)  be  linear  in  Reynolds  stress;  (ii)  confirm  to  a  specific 
functional  form  (given  in  equation  1 .2.6);  and,  (iii)  lead  to  realizable  Reynolds  stress.  While  the 
need  for  the  realizability  requirement  in  now  clear,  the  linearity  requirement  also  stems  from 


important  turbulence  physics.  It  is  easy  to  see  that  Poisson  equation  for  rapid  pressure  is  linear  in 
fluctuating  velocity  (Reynolds,  1976).  Therefore,  the  rapid  or  linear  pressure  term  cannot 
increase  or  proliferate  the  number  of  (Fourier)  modes  of  fluctuation.  In  other  words,  the  rapid 
pressure  term  conserves  the  number  of  fluctuating  modes.  Any  non-linear  model  cannot  preserve 
this  important  physical  character  rapid  pressure.  The  Miju-form  requirement  is  a  consequence  of 
the  form  of  the  Greens  function  solution  of  the  Poisson  equation.  As  pointed  out  in  Pope  (2000), 
currently,  there  exists  no  rapid  pressure-strain  correlation  model  that  satisfies  all  the  required 
constraints.  Specifically,  none  of  the  current  models  satisfy  all  realizability  and  linearity 
requirements.  In  fact,  there  has  been  much  speculation  that  it  is  virtually  impossible  to  satisfy 
realizability  fully  wifli  linear  models  (Lumley,  1978).  While  realizability  is  an  important 
constraint,  it  is  equally  desirable  to  achieve  it  in  a  manner  consistent  with  turbulent  physics,  that 
is,  without  compromising  on  the  other  two  physical  requirements. 


In  this  paper  we  propose  additional  constraints  on  0'.  and  Myiu  that  should  enable  the 
development  of  more  accurate  closure  models. 


1.3  Realizability  constraints 

We  will  first  discuss  the  strong  and  weak  forms  of  Schumann  constraints.  Realizability  of 
Reynolds  stresses  is  most  conveniently  considered  in  terms  of  the  determinant  of  the  normalized 
Reynolds-stress  tensor: 


(1.3.1) 


In  the  realizable  state-space  of  Reynolds  stress  (Lumley  invariant  triangle),  Fis positive:  F  is 
negative  through  most  of  the  unrealizable  space.  The  realizable  and  the  unrealizable  Reynolds 
stress  state-spaces  are  separated  by  two-componential  turbulence  in  which  determinant  F 
vanishes.  The  current  realizability  conditions  constrain  die  behavior  of  the  model  in  the 
neighborhood  of  F=0. 


Strong  realizability.  The  strong  Reynolds  stress  realizability  condition  requires  the  turbulence 
model  to  satisfy  the  following  conditions  (Pope,  2000): 

,dF,  .  /F.  . 

(^>■"=0;  (^>-'»o.  (1.3.2) 

The  strong  realizability  condition  permits  accessibility  to  two-component  turbulence.  Rigorous 
implementation  of  strong  realizability  constraints  to  develop  closure  models  is  best  exemplified 
in  Ristorcelli,  Lumley  and  Abid  (1995),  Sjogren  and  Johansson  (2000),  etc.  In  fact,  in  Johansson 
and  Hallback  (1994)  a  recipe  for  deriving  a  rapid  pressure-strain  correlation  model  that  satisfies 
the  strong  realizability  constraint  at  any  order  of  tensor  expansion  is  provided.  It  is  found  that  the 
fourth-order  model  was  successful  over  a  wider  range  of  flows  that  typical  linear  models.  In 
Sjogren  and  Johansson  (2000),  the  use  of  strong  realizability  constraint  to  develop  slow  pressure 
correlation  and  anisotropic  dissipation  are  also  siunmarized. 


Weak  realizabUity.  Many  authors  point  to  major  difficulties  in  enforcing  the  strong  realizability 
constraint  and  opt  for  the  easier  to  implement  weak  realizability  version. 


The  weak  reali2ability  constraint  simply  requires  that  the  rate  of  change  of  F  be  positive  as  the 
two-componential  limit  is  approached  (Pope,  2000): 

( - )f  =  o>0.  (1.3.3) 

dT 

This  constraint  doesn’t  permit  access  to  two -component  turbulence.  The  stochastic  realizability 
analysis  of  Durbin  and  Speziale  (1994)  is  a  variant  of  the  weak-realizability  approach.  They 
demonstrate  that  every  turbulence  model  equation  can  be  represented  by  an  equivalent  stochastic 
equation.  The  requirement  that  the  coefficient  of  the  random  forcing  term  be  real  yields 
constraints  on  the  pressure-strain  conelation  model.  This  approach  leads  to  the  dependence  of 
the  rapid  pressure-strain  correlation  model  coefficients  on  the  coefficients  in  the  dissipation  and 
the  slow  pressure- strain  c  orrelation  models.  Such  interdependence  of  model  coefficients  is  not 
physically  justifiable.  F  or  example,  if  turbulence  model  yields  unrealizable  results  in  the  rapid 
distortion  limit,  it  is  the  rapid  pressure-strain  correlation  model  that  needs  modification.  Slow 
pressure- strain  correlation  and  dissipation  processes  are  completely  irrelevant  in  fiiis  limit  and 
modifying  their  models  to  achieve  realizability  must  be  deemed  unphysical.  The  realizability 
issues  of  each  process  should  be  considered  independently  as  is  well  demonstrated  in  Sjogren 
and  Johansson  (2000). 

Need  for  further  constraints.  In  the  second-moment  closure  approach,  Reynolds  stress  is 
calculated  using  a  modeled  evolution  equation.  Any  non-physical  behavior  of  fire  Re)molds 
stress  must  originate  from  unphysical  closure  models  in  its  evolution  equation.  In  Navier-Stokes 
turbulence,  the  realizability  of  the  Reynolds  stress  is  achieved  while  maintaining  all  the 
contributing  statistical  moments  individually  realizable.  Any  realizability  analysis  is  incomplete 
without  the  examination  of  the  realizability  of  the  individual  dosure  models  that  contribute 
towards  Reynolds  stress  evolution.  The  strong  realizability  approach  indirectly  addresses  this 
issue,  but  only  in  the  two-componential  turbulence  limit.  The  weak  realizability  approach,  as 
currently  practiced,  completely  ignores  the  realizability  of  individual  moments  in  the  Reynolds 
stress  evolution  equation.  Thus,  both  the  strong  and  weak  realizability  approaches  can  and  do  (as 
will  be  demonstrated  later)  permit  unrealizable  pressure-strain  correlation  values. 

In  order  to  achieve  a  greater  degree  of  consistency  with  turbulence  physics,  we  recommend  that 
the  closure  model  for  each  moment  in  the  Reynolds  stress  evolution  equation  be  individually 
realizable,  not  only  in  the  two-componential  limit  but  also  at  all  other  states  of  turbulence. 
Physically  achievable  or  realizable  models  for  each  of  the  unclosed  statistical  moments  must 
satisfy  the  following  requirements: 

(a)  Auto-covariances  must  be  non-negative. 

(b)  Cross-covariances  must  satisfy  Schwarz  inequality. 

In  this  paper,  we  address  the  realizability  of  the  rapid  pressure-strain  correlation  term.  In  most 
modeling  mefliods,  it  is  the  Myki  tensor  that  is  first  modeled  -  see  for  example.  Launder,  Reece 
and  Rodi  (1975);  Ristorcelii,  Lumley  and  Abid  (1995)  or  equation  (54)  in  Sjogren  and 
Johansson  (2000)-  and  file  pressure- strain  correlation  closure  is  simply  obtained  from  equation 
(1.2.6).  Hence,  postulating  constraints  on  Afy*/  is  as  effective  and  useful  as  constraining 
pressure-strain  correlation  itself  Here,  we  derive  constraints  on  the  Miju  tensor  that  guarantees 
realizable  rapid  pressure-stress  correlation. 


The  need  for  further  constraints  on  Myia  can  also  be  motivated  from  the  turbulence  structure 
tensors  point  of  view.  As  Kassinos  et  al  (2001)  point  out,  a  statistically  complete  one-point 
description  of  turbulence  requires  specification  several  tensors:  componentiality  or  Reynolds 
stress  (R,y);  Dimensionality  (Dj,);  Circularity  (Cj,);  and,  finally,  stropholysis  (0;/*).  In  general, 
each  tensor  carries  independent  information  and  all  of  them  are  needed  to  fully  specify  the  one- 
point  statistical  state  of  turbulence.  It  then  stands  to  reason  that  closure  for  all  tensors  should  be 
individually  realizable.  In  the  special  case  of  incompressible  homogeneous  turbulence  these 
structure  tensors  can  be  related  to  the  Myki  tensor  (Kassinos  et  al,  2001): 

Qijk  ~  ^ipq  ^jqpkt  Rij  ~  ^ijrr  /  Dy  Mrry  Fy—  Bimp  Bjrs^psrm  >  (I  •3*4) 

where  eyk  is  the  alternating  tensor.  Therefore,  in  homogeneous  incompressible  turbulence,  all 
physical  and  realizability  issues  are  embodied  in  Myu  tensor.  Hence,  constraints  that  can 
guarantee  the  physical  fidelity  of  this  tensor  are  of  great  importance. 

1.4  Realizability  of  rapid  pressure  moments 

The  concept  of  realizability  requires  that  the  closure  models  for  turbulence  statistics  are  based  on 
a  velocity  field  that  is  realizable  or  physically  achievable.  But  what  is  a  physically  possible 
fluctuating  velocity  field?  It  is  generally,  although  tacitly,  assumed  that  any  real- valued  vector 
field  (as  opposed  to  complex  or  imaginary  fields)  is  an  achievable  velocity  field.  Any  real- valued 
vector  field  would  yield  non- negative  values  for  the  diagonal  components  of  Reynolds  stress. 
Only  complex  or  imaginary  velocity  field  would  lead  to  negative  energy  components.  This 
argument  is  the  basis  of  Schumaim  constraints  which  permit  real  vector  fields. 

We  argue  that  not  all  real- valued  vector  fields  are  achievable  velocity  fields.  Additional 
kinematic  constraints  must  be  satisfied  including  the  compressibility  condition.  If  a  flow  is 
incompressible,  then  the  velocity  vector  field  must  be  divergence  free  at  all  points  in  space  and 
time.  Equivalently,  the  pressure  field  obtained  by  Poisson  equation  from  a  candidate  real- valued 
vector  field  must  also  be  real- valued.  Thus,  a  fully  realizable  turbulent  velocity  field  must  not 
only  be  real-valued,  but  also  be  divergence  free.  Alternatively,  the  implied  pressure  field  must  be 
real.  The  incompressibility  aspect  of  the  velocity  field  is  not  considered  in  any  of  the  current 
realizability  methods.  While  fire  Schumarm  constraints  are  necessary  for  achieving  a  realizable 
turbulence  model,  they  are  not  sufficient. 

The  incompressibility  requirement  can  be  enforced  by  constraining  the  pressure  fluctuations.  It 
id  the  pressure  field,  after  all,  tiiat  preserves  the  incompressibility  of  the  flow.  Further,  pressure 
gradient  is  responsible  of  the  evolution  of  the  velocity  field.  It  is,  therefore,  essential  to  mandate 
realizable  pressure  statistics  along  with  constraints  on  the  Reynolds  stress.  Ibr  compressible 
velocity  fields,  the  divergence  of  the  velocity  field  must  be  consistent  with  density  fluctuations. 
Then,  the  realizability  issues  of  the  velocity  field  will  be  inextricably  coupled  with  those  of 
pressure  and  temperature  statistical  moments.  Compressible  flows  will  be  considered  in  future 
works. 

Realizability  bounds.  In  homogeneous  turbulence,  the  pressure- strain  correlation  (^^.)  is 
identical  in  magnitude  and  opposite  in  sign  to  the  velocity  pressure- gradient  correlation  ( n,^): 


(1.4.1) 


^ij  =  ^  =-2Ay  =-0;- 
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dx, 


The  velocity  pressure- gradient  correlation  is  more  fundamental  as  it  is  the  one  that  appears  in  the 
Reynolds  stress  closure  equation  before  decomposition  into  homogeneous  and  inhomogeneous 
components.  We  will  now  investigate  the  realizability  bounds  on  the  rapid  portion  of  the  velocity 
pressure-strain  correlation  (  \j).  The  magnitude  of  this  correlation  is  bound  by 
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(1.4.2) 


Schwarz  inequality  as  applied  to  the  last  two  correlations  on  the  extreme  right  requires 

Sp'\  ' 
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(1.4.3) 


(Whereas  repeated  English  indices  imply  summation,  Greek  indices  do  not.)  As  a  consequence, 
the  rapid  portion  of  velocity  pressure  correlatfan  must  be  bounded  by; 
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Since— ,  the  following  inequality  must  hold: 
^  ^  dx,  dx. 
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where  the  second  identity  defines  the  tensor  .  Due  to  the  fact  that  we  replace  — r~  in  the 
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inequality  by  the  larger-^ — —  ,  we  must  regard  the  constraint  in  equation  $.4.5)  as  less 

dx,  dx, 

stringent  than  required  to  satisfy  the  Schwarz  inequality.  This  inequality  is  valid  for 
homogeneous  and  inhomogeneous  flows.  From  equation  (1.4.5),  it  can  be  seen  that  the 
realizability  bounds  on  the  rapid  pressure-strain  correlation  depends  on  the  componentiality  (via 

_ 2. 

-  dp'  2 

u„u„ )  and  dimensionality  (via - )  of  the  turbulence. 

dxp  dxp 


The  inequality  in  expression  0  -4.5)  forms  the  basis  of  the  new  constraints  developed  in  this 
paper.  In  the  limit  of  two-componential  turbulence,  this  inequality  yields  the  strong  realizability 
constraint:  if  velocity  fluctuations  in  the  a-direction  vanish,  then 


=0, 

which  immediately  implies 
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Apart  from  this  requirement,  the  inequality  in  equation  0-4.5)  has  not  been  exploited  for 
formulating  all  other  possible  modeling  constraints.  For  example,  it  is  obvious  that  when  the 
pressure- gradient  variance  vanishes,  so  must  die  pressure-strain  correlation.  Yet,  this 
fundamental  requirement  has  never  been  posed  as  a  modeling  constraint.  This  is  most  likely  due 
to  the  fact  that  rapid  pressure- gradient  statistics  (indicative  of  the  dimensionality  of  turbulence) 
are  not  considered  in  traditional  modeling  approaches  where  only  the  Reynolds  stresses 
(componetiality  of  turbulence)  is  solved  for  explicitly.  In  this  work,  we  will  derive  an  exact 
analytical  expression  for  rapid  pressure-gradient  variance.  Knowledge  of  rapid  pressure-gradient 
variance  opens  up  the  possibility  of  a  new  set  of  constraints  of  the  rapid  pressure-strain 
correlation  that  are  much  more  general  that  the  strong  realizability  requirement. 


Rapid  pressure -gradient  variance.  Multiplying  equation  (1.2.3)  through  by  the  rapid  pressure 
and  taking  average  leads  to: 


In  homogeneous  turbulence  which  is  of  interest  here,  we  can  further  write 
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It  was  shown  in  the  previous  section  that  the  rapid  pressure-strain  correlation  tensor  is  related  to 
the  fourth  order  tensor  Mjki.  The  relationship  between  this  tensor  and  the  pressure-gradient 
variance  can  be  easily  derived  from  equations  (1 .2.4)  and  (1 .4.9): 
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Thus,  the  rapid  pressure-gradient  variance  can  be  obtained  as  a  function  of  the  rapid  pressure- 
strain  correlation.  In  purely-strained  flows,  the  relationship  is  even  simpler: 


p’-^-^  =  2S,,p's„ 

^  dxdx 


(1.4.11) 


It  has  to  be  pointed  out  that  the  Kassinos,  Reynolds  and  Rogers  (2001)  derive  an  expression  for 
pressure- gradient  variance  in  die  rapid  reference- frame  rotation  limit.  The  expression  involves 


tile  rotation  rate  and  the  circularity  tensor  which  contains  information  about  the  large-scale 
vorticity.  The  circularity  tensor  is  not  known  the  traditional  second-moment  closure  modeling. 

New  realizability  constraints 


With  the  knowledge  of  rapid  pressure-gradient  variance,  we  will  now  formulate  the  new 
realizability  constraints  in  one-,  two-  and  three-componential  homogeneous  turbulence.  These 
pressure-strain  correlation  constraints  account  for  dimensionality  as  well  as  componetiality, 
whereas,  the  Schumann  constraints  takes  only  componentailty  into  consideration. 


Constraint  #1.  In  two-componential  turbulence  (with  no  velocity  fluctuations  in  the  a  direction), 
the  realizability  of  rapid  pressure- strain  correlation  requires 

The  values  of  tiie  off-diagonal  terms  will  depend  on  the  dimensionality  of  the  turbulence.  This  is 
similar  to  the  strong  realizability  requirement.  Alternatively,  the  realizability  constraint  on  Myu 
in  the  two-component  limit  is 


Maaij=0.  (1.4.13) 

In  one-component  turbulence,  all  three  diagonal  pressure-strain  correlation  components  must  be 
zero.  The  diagonal  components  in  the  direction  of  zero  fluctuating  velocity  clearly  must  be  zero 
for  the  same  reason  as  above.  The  correlation  component  in  the  third  direction  is  zero  for  a 
different  reason.  Continuity  requires  that  the  fluctuating  pressure  gradient  be  zero  in  the  direction 
of  velocity.  Thus,  in  one-dimensional  turbulence  (if  6  is  the  direction  of  non- zero  velocity)  the 
realizability  constraints  are 


Vox  =  0,  implying  =  0;  and  Maaij  =  0. 

For  a  =  =  0 ,  implying  <p'^  =  0;  and  Mijaa =0.  (1 .4. 14) 

d)c„  dx„ 

The  pressure-strain  constraints  are  valid  for  the  slow  term  as  well. 


Constraint  #  2.  The  rapid  pressure-strain  correlation  model  must  lead  to  non-negative  values  for 
the  rapid  pressure- gradient  variance.  Physically,  this  ensures  that  the  pressure  fluctuations  are 
real  and,  hence,  the  velocity  field  is  divergence- free.  Thus,  the  first  of  two  realizability 
constraints  on  Myid  is 


a.,  a,  ■”*- 

This  inequality  should  hold  for  any  arbitrary  non-zero  mean- velocity  gradient  tensor.  Hence,  the 
Mjjki  must  be  positive  semi-definite.  Unless  Mykt  is  positive  semi-definite,  tiiere  will  be  a  mean 
velocity  field,  for  which  the  fluctuating  rapid  pressure  gradient  will  be  negative. 


For  purely  framed  flows,  positivity  condition  on  the  pressure-gradient  variance  can  also  be 
written  as 


(1.4.16) 


Constraint  #  3  The  pressure-gradient  variance  must  not  only  be  positive,  but  also  be  large 
enough  to  satisfy  tiie  Schwarz  inequality  on  rapid  pressure-strain  correlation.  This  can  be 
guaranteed  if  the  model  satisfies  the  following  inequality  (firom  equations  1.4.10  and  1.4.5): 
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Satisfaction  of  tiiis  constraint  in  conjunction  with  constraint  #  2  will  ensure  realizable  velocity 
pressure- gradient  correlation.  For  purely  strained  flows  the  above  constraint  can  be  expressed 
directly  in  terms  of  pressure-strain  correlation: 
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(1.4.18) 


It  can  be  recognized  that  the  first  constraint  (equation  1.4.12),  derived  for  one-and  two- 
componential  turbulence,  is  equivalent  to  the  strong  realizability  constraint.  However,  constraint 
#  2  (equation  1.4.15)  and  #3  (equation  1 .4.17)  are  not  addressed  in  existing  literature.  In  light  of 
these  new  constraints  several  questions  arise: 

(a)  Do  the  current  models  violate  the  new  realizability  constraints?  If  not,  the  new  constraints 
are  only  of  academic  interest. 

(b)  Does  Reynolds-stress  realizability  automatically  imply  realizability  of  pressure- strain 
correlation?  If  yes,  then  tfie  current  constraints  would  be  redundant. 

(c)  Conversely,  does  pressure-strain  correlation  realizability  guarantee  Reynolds-stress 
realizability?  If  yes,  then  the  new  realizability  constraints  are  more  complete  and  the 
Reynolds-stress  realizability  is  redundant. 

(d)  Or,  are  the  old  and  new  realizability  constraints  mutually  exclusive  so  that  both  should  be 
required  of  the  pressure-strain  correlation  models? 

In  an  attempt  to  find  the  answers,  we  perform  model  calculations  in  the  next  section. 


1.5  Solution  procedure 

Fourth-order  Runge-Kutta  scheme  is  used  for  time  integration  of  second-order  closure  equations. 
The  mean  flows  investigated  are 
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Evolution  trajectories  are  tracked  on  the  invariant  plane,  where  the  invariants  are  defined 
by  6t)^  =bybj,  and  6^^=bgb^bu.  For  each  mean  flow,  numerous  initial  values  of  {^.r}) 

encompassing  the  entire  Lumley  triangle,  are  chosen.  Each  point  on  flie  Lumley  triangle 
corresponds  to  six  combinations  of  anisotropy  tensor  components.  All  six  combinations  are 
investigated.  All  together,  it  results  h  over  64000  initial  conditions  considered  for  each  mean 
flow.  In  all  simulations,  initial  values  of  k  (ko)  and  e  (£<>)  are  taken  (without  any  loss  of 

generality)  to  be  unity.  Calculations  are  performed  for  three  values  of  S*  =  Skje^ :  5,  15  and  50. 

The  case  S*  =5  corresponds  to  near  equilibrium  turbulence,  while  5*  =50  corresponds  to  rapidly 
distorted  turbulence. 

1.6  Realizability  study 

The  results  are  presented  in  four  steps:  (i)  the  set  of  anisotropy  initial  conditions  that  ultimately 
lead  to  RS-realizability  violation  are  identified  for  each  model;  (ii)  contribution  of  the  various 
terms  in  the  Reynolds  stress  evolution  equation  towards  unrealizable  trajectories  is  examined, 
and  the  term  responsible  for  RS-realizability  violation  is  identified;  (iii)  connection  between  RS- 
realizability  violation  and  an  unrealizable  rapid  pressure- strain  correlation  model  is  analyzed; 
(iv)  tmncated  LRR  model  (Girimaji,  2002),  which  satisfies  all  constraints  listed  in  Pope  (2000) 
and  also  rapid  PSC-  realizability  constraint  (1.4.17),  is  considered. 

1.6.1  RS-realizability  violation 

Solution  trajectories  obtained  firom  model  calculations  can  be  classified  into  two  types,  namely 
the  ones  that  are  always  realizable  and  asymptote  to  equilibrium  and  the  ones  that  violate 
realizability  (leave  the  Lumley  triangle).  The  area  in  the  Lumley  triangle  containing  initial 
conditions  fi-om  which  only  realizable  trajectories  originate  is  called  the  basin-of-attraction.  The 
area  inside  the  Lumley  triangle,  but  outside  the  basin-of-attraction  will  be  called  as  the  basin-of- 
violation.  Notice  though,  that  initial  conditions  inside  the  basin-of-violation  may  or  may 
not  lead  to  uiu’ealizable  solution  trajectories  depending  on  the  combination  of  b^  -components. 
The  line  that  demarcates  two  basins  is  called  the  violation  boundary.  The  location  of  the 

violation  boundary  depends  upon  pressure-strain  model  employed,  mean  flow,  and  5*  -value.  All 
initial  conditions  below  the  violation  boundary  belong  to  the  basin-of-attraction. 

Violation  boundaries  for  the  LRR-QI  model  are  shown  in  figure  1  at  S*  -S  and  in  figure  2  at 

5*  =  50 .  At  5*  =  5 ,  the  smallest  basirr-of- violation  is  observed  in  the  homogenous  shear  flow  and 
the  largest  one  in  the  axisymmetric  expansion  flow.  Basins-of-violation  for  other  three  cases  fall 
between  those  two  regions.  The  basinrof-violation  in  the  axisymmetric  contraction  flow  is 
smaller  than  the  one  in  the  axisymmetric  expansion  flow;  and  such  region  in  the  plain  strain  flow 
is  almost  similar  to  the  one  in  the  mixed  shear  flow.  The  mixed  shear  flow  considered  here  has 
80%  of  strain  and  only  20%  of  rotation.  Hence,  basins-of-violation  in  this  flow  and  in  the  plain 
strain  flow  are  similar.  Once  the  ratio  of  strain  to  rotation  is  decreased,  the  mixed  shear  basins- 
of  violation  start  to  transform  towards  the  basins-of-violation  corresponding  to  the  homogeneous 

shear  flow.  With  increasing  S* ,  a  violation  boundary  goes  down,  and  the  basin-of-violation  area 
increases.  The  overall  shape  of  the  region  does  not  change  though.  Such  tendencies  are  observed 
in  all  flows  (see  figures  1.1  and  1.2). 


The  unrealizable  regions  corresponding  to  the  SSG  model  (figures  1.3  and  1.4)  behave  similarly. 
The  only  notable  difference  between  LRR-QI  and  SSG  models  is  that  basins-of- violation  for  the 
SSG  model  are  larger  than  those  of  die  LRR-QI  model  in  all  flows  considered. 

Selected  solution  trajectories  evolving  from  the  basin-of-attraction  and  outside  that  area  are 
presented  in  figures  1.5- 1.9.  In  the  figures,  a  violation  boundary  is  shown  by  solid  thick  line.  The 
cross  (@)  marks  the  equilibrium  point  in  the  flow.  The  values  of  the  anisotropy  tensor 
components  in  the  equilibrium  point  are  shown  in  table  1.1  for  different  flow  conditions  and  for 
two  models:  LRR-QI  and  SSG.  Some  of  equilibrium  -values  given  in  table  1.1  slightly  differ 

from  those  used  in  previous  studies.  It  is  explained  by  different  C^j  value  employed  in  those 
studies.  The  rest  of  the  symbols  in  figures  1.5- 1.9  show  initial  conditions  from  which  solution 
trajectories  evolve.  Arrows  on  trajectories  show  the  direction  of  trajectory  evolutioa 


Table  1.1  Equilibrium  anisotropy  values. 


LRR-QI  model 

bii  b22 

b33 

bi2 

SSG  model 

b„ 

1)22 

b33 

bi2 

A 

C 

0.2513 

0.1256 

0.1256 

0.0000 

0.2267 

0.1134 

0.1134 

0.0000 

HS 

0.1565 

0.1223 

0.0342 

0.1877 

0.2216 

0.1485 

0.0731 

0.1640 

AE 

0.2855 

0.1428 

0.1428 

0.0000 

0.3216 

0.1608 

0.1608 

0.0000 

PS 

0.2167 

0.2510 

0.0342 

0.0000 

0.1861 

0.2548 

0.0688 

0.0000 

M 

S 

0.1992 

0.2334 

0.0342 

0.0889 

0.1885 

0.2616 

0.0731 

0.1025 

In  the  homogenous  shear  flow  (figure  1.5)  trajectories  evolving  from  initial  conditions  below  the 
violation  boundary  end  up  in  the  equilibrium  point.  In  the  plain  strain  flow,  realizable  trajectories 
evolve  in  the  complicated  manner  (figure  1 .6)  as  they  make  loops  before  reaching  equilibrium. 
The  behavior  of  the  realizable  trajectories  in  the  mixed  shear  flow  (figure  1.7)  resembles  the 
behavior  of  trajectories  in  both  plain  strain  and  homogenous  shear  flows.  In  the  presence  of  shear 
though,  loops  are  reduced  and  they  do  not  reach  axisymmetric  expansion  and  axisymmetric 
contraction  limits,  as  it  would  be  in  the  plain  strain  flow  (figure  1 .6).  The  shape  of  trajectories 
appears  to  depend  on  the  ratio  of  strain  to  shear.  It  is  observed  that  when  the  ratio  reduces, 
trajectories  evolve  closer  to  the  homogenous  shear  flow  case.  In  the  axisymmetric  contraction 
flow,  die  trajectories  (?,  |  ,♦)  that  asymptote  to  equilibrium  also  exhibit  the  looping  behavior 
(figure  1.8).  In  the  axisymmetric  expansion  flow,  similar  behavior  of  realizable  trajectories  is 
observed  (figure  1.9).  In  figures  1.5- 1.9,  trajectories  which  start  from  initial  conditions  above  the 
violation  boundary  evolve  outside  the  Lumley  triangle  and  so,  are  unrealizable. 

Results  of  this  subsection  show  that  there  exists  a  sizable  finite  set  of  initial  conditions  for  which 
the  solution  trajectories  violate  realizability.  The  basin-of-violation  grows  with  increasing  the 


parameter  S*  in  all  flows.  It  also  depends  on  the  rapid  pressure-strain  model  and  the  mean  flow. 
We  will  now  tiy  to  identify  which  term  in  the  Reynolds  stress  evolution  equation  contributes 
most  toward  RS -realizability  violation. 

1.6.2  Closure  model  responsible  for  RS-realizability  violation 

The  closure  model  responsible  for  the  RS-realizability  violation  can  be  identified  by  examining 
the  contribution  of  various  terms  in  the  Reynolds  stress  evolution  equation.  The  evolution 
equation  of  the  Reynolds  stresses  can  be  split  into  three  parts: 


1.  due  to  production 


2.  due  to  rapid  pressure-strain 


\d  — 
Idt  '  I 


3.  due  to  non-linear  processes  of  dissipation  and  slow  pressure-strain  correlation 


(1.6.1) 


(1.6.2) 


At  every  point,  the  resultant  firom  these  three  parts  will  determine  the  direction  of  a  trajectory. 
Some  umealizable  trajectories  are  shown  in  figures  1.10  and  1.11  in  plain  strain  and 
axisymmetric  expansion  flows.  They  correspond  to  the  initial  conditions  (a)  in  figure  1 .6  (PS) 
and  (*)  in  figure  1 .9  (AE).  Contribution  of  the  production,  the  rapid  pressure-strain  correlation, 
and  non-linear  processes  in  trajectory  evolution  is  depicted  at  two  different  stages,  when  Miju  and 
RS  violations  occur. 

In  figure  1.10,  contribution  of  different  parts  towards  anisotropy  evolution  is  qualitatively  similar 
at  both  points.  Details  can  be  summarized  as  follows: 

•  Production  ^)  drives  the  trajectory  parallel  to  the  two-component  limit  and  hence,  does 
not  lead  to  realizability  violation.  Since  this  term  is  ‘exact’,  it  is  indeed  expected  to  be 
physically  realizable. 

•  Slow  pressure- strain  and  dissipation  (s)  always  attract  the  trajectory  towards  the  interior 
of  the  Lumley  triangle.  These  terms  have  isotropizing,  rather  than  anisotropizing  effect. 

•  Rapid  pressure- strain  correlation  (r)  points  in  a  direction  away  from  the  Lumley  triangle. 
It  is  clear  that  this  is  the  component  responsible  for  RS  realizability  violation. 

The  same  trends  are  seen  in  figure  1.11  and  in  other  umealizable  trajectories  investigated 
(figures  not  shown). 

It  can  be  said  with  a  great  degree  of  generality  that  this  is  a  rapid  pressure-strain  correlation 
model,  which  leads  to  an  umealizable  trajectory  of  the  Reynolds  stresses  in  considered  flows. 
Thus,  conclusion  of  this  subsection  is  that  the  rapid  pressure-strain  component  is  responsible  for 
Reynolds  stress  realizability  violation. 


1.6.3  Connection  between  RS-  and  rapid  PSC-reaiizabilities 

All  of  die  unrealizable  trajectories  from  the  previous  subsection  are  now  investigated  for 
compliance  with  Myu  constraint  (1.4.17).  As  it  is  seen  on  figures  1.10  and  1.11,  violation  of  the 
Mijki  constraint  precedes  RS -realizability  violation.  It  happens  in  all  cases  considered.  In  figures 
1.12-1.15,  the  line  formed  by  points  at  which  theM,yt/  constraint  is  first  violated  by  unrealizable 
trajectories,  is  shown  in  various  mean  flows  for  LRR-QI  and  SSG  rapid  pressure- strain 
correlation  models.  We  will  call  this  line  as  the  violation  boundary.  In  all  these  cases  the 
line  is  nearly  parallel  to  the  two-component  line  of  the  Lumley  triangle.  Comparing  SSG  and 
LRR-QI  models,  one  can  say  that  the  former  violates  the  Mya  constraint  over  larger  portion  of 
the  Lumley  triangle  dran  the  LRR-QI  model. 

Investigation  of  realizable  trajectories  reveals  a  surprising  fact.  Although  most  of  them  satisfy 
the  Miju  constraint,  some  trajectories  do  violate  it.  In  figure  1.16,  two  sample  trajectories 

calculated  with  the  LRR-QI  model  are  shown  in  the  axisymmetric  expansion  flow  at  5*  =50 . 
One  of  the  trajectories  (starting  from  ?)  becomes  umealizable  at  (®)  after  violating  the  Mya 
constraint  at  (^).  The  other  trajectory  (starting  from  s)  violates  the  Myki  constraint  at  (^),  but  the 
Reynolds  stresses  remain  realizable  at  all  times.  In  fact,  at  later  stages  dre  Myu  constraint  is  again 
satisfied  by  the  very  same  trajectory  as  it  asymptotes  to  the  equilibrium  point.  However,  because 
the  Af{,A/-constraint  violation  occurs,  the  behavior  of  this  trajectory  must  be  deemed  unphysical. 

Figure  1.16  demonstrates  clearly  that  realizable  Reynolds  stresses  do  not  always  guarantee  the 
realizable  rapid  pressure-strain  correlation.  It  means  that  new  Mya  constraint  (1.4.17)  carmot  be 
guaranteed  even  by  the  strong  realizability  condition.  Prompted  by  this  finding,  an  investigation 
of  noiylinear  rapid  pressure-strain  correlation  models  (Shih  and  Lumley,  1985;  Sjogren  and 
Johansson,  2000)  that  satisfy  the  strong  realizability  condition  is  performed.  The  SL  model 
always  produces  realizable  trajectories.  However,  this  model  yields  unrealizable  rapid  pressure- 
strain  correlations  over  a  significant  portion  of  the  Lumley  triangle  near  the  two-component  line. 
These  uiurealizable  zones  for  the  plain  strain  and  axisymmetric  contraction  cases  are  shown  in 
figure  1.17. 

Calculations  were  also  performed  with  a  fourth-order  rapid  pressure-strain  correlation  model  of 
Sjogren  and  Johansson  (2000).  The  coefficients  of  this  model  are 

g,  = -0.1 43,gj  =00295,^3  =-00484,g,  =14. 

The  relationship  between  g’s  and  f  s  can  be  found  in  Sjogren  and  Johansson  (2000).  It  was  found 
that  this  model  yields  fully  realizable  results  everywhere  except  the  near  vicinity  of  the  two- 
component  turbulence  line.  While  the  Reynolds  stresses  are  always  realizable  (as  guaranteed  by 
the  model  development),  tiie  pressure-strain  correlation  shows  small  region  of  imrealizability. 

The  final  conclusion  from  this  subsection  is  that  the  rapid  pressure-strain  realizability  constraint 
cannot  be  guaranteed  by  compliance  with  strong  realizability. 

1.6.4  Truncated  LRR  model 

The  new  PSC-realizability  constraint  ^.4.17)  should  aid  in  the  development  of  simple  rapid 
pressure-strain  correlation  models  with  higher  degree  of  fidelity  to  turbulence  physics.  As  a 


simple  demonstration,  a  new  truncated  (clipped)  LRR  model  for  the  pressme-strain  correlation 
has  been  derived  in  Girimaji  (2002): 


[R^j  +  R^j  if  Constraint  (1.4.17)  is  valid 
\Ry  if  (1 .4. 17) is  violated 


(1.6.4) 


In  the  model,  terms  R^j  and  Rij  are  corresponding  slow  and  rapid  components  of  the  LRR-QI 
pressure- strain  correlation  model.  This  model  is  piece-wise  linear.  Due  to  its  simplicity 
(compared  to  non-linear  models  (Shih  and  Lumley,  1985;  Sjogren  and  Johansson,  2000) 
satisfying  strong  realizability),  it  is  likely  to  be  of  hi^  practical  value.  We  do  not  advocate  the 
use  of  fliis  clipped  model  in  practical  computations.  Rather,  we  use  it  here  to  demonstrate  that 
satisfying  the  new  constraints  will  guarantee  Reynolds  stress  realizability. 

By  testing  model  (1.6.4)  in  different  flows,  it  was  confirmed  that  the  truncated  LRR  model  is 
realizable  for  all  initial  conditions  in  file  Lumley  triangle.  Sample  corrected  trajectories  are 

shown  in  the  plain  strain  flow  at  S*  =5  (figure  1.18)  and  in  the  axisymmetric  expansion  flow  at 

S*  =  15  (figure  1.19).  In  figure  1.18  (a),  the  trajectory  obtained  with  the  original  LRR-QI  model 
(solid  line)  starts  from  the  initial  condition  Q).  Then,  it  violates  the  rapid  PSC-realizability 
constraint  at  (^)  and  fire  RS -realizability  at  Q).  This  trajectory  is  corrected  in  the  truncated 
model  (dash  line).  It  is  clearly  seen  that  the  trajectory  returns  to  equilibrium  after  taking  a 
circuitous  and  fully  realizable  path.  The  same  is  observed  on  figures  1.18  (b)  and  1.19  (a,b).  In 
Figures  1.20  and  1.21  we  show  various  contributions  toward  the  evolution  of  original  and 
truncated  LRR  model  trajectories.  The  contributions  are  shown  at  two  different  stages:  RS-  and 
rapid  PSC-realizability  violation  states.  Initial  conditions  are  same  as  Arose  in  figures  1.18  (b) 
and  1.19(b). 

1.7  Summary  and  Conclusions 

A  detailed  study  of  Reynolds  stress  realizability  issues  is  presented  in  this  paper.  It  is  shown  that 
all  realizable  trajectories  asymptote  to  a  unique  equilibrium  point.  The  set  of  anisotropy  initial 
conditions  that  ultimately  1  ead  to  Reynolds  stress  realizability  violation  (basin-of- violation)  is 
identified  for  each  model  (Speziale  et  al,  1991;  Launder  et  al,  1975;  Shih  et  al,  1985)  in 

different  mean  flows  and  at  different  initial  values  of  5*.  The  basin-of- violation  grows  with 

increasing  the  parameter  S* .  It  also  depends  on  the  rapid  pressure-strain  correlation  model  and 
the  mean  flow  considered. 

Contribution  of  various  terms  of  the  model  Reynolds  stress  equation  toward  trajectory  evolution 
is  studied.  The  rapid  part  of  the  pressure-strain  correlation  is  identified  to  be  the  term  responsible 
for  the  Reynolds  stress  realizability  violation.  In  this  connection,  the  realizability  of  a  rapid 
pressure- strain  correlation  closure  is  examined,  and  its  direct  coimection  with  RS -realizability  is 
established  in  the  first  time.  Computations  demonstrate  that  violation  of  the  rapid  pressure-strain 
correlation  realizability  precedes  violation  of  Rejmolds  stress  realizability  in  all  considered  cases. 
Also,  Reynolds  stress  realizability  based  on  either  weak  or  strong  formulation  does  not  guarantee 
the  realizable  rapid  pressure- strain  correlation.  Thus,  realizability  constraint  (1.4.17)  on  flie  rapid 


pressure-strain  correlation  is  more  fundamental.  Such  understanding  of  the  turbulence  model 
dynamics  is  cracial  for  developing  a  realizable  and  physically  consistent  turbulence  model. 

Using  die  rapid  pressure-strain  correlation  realizability  constraint  (1.4.17),  the  piece- wise  linear 
model  (1.6.4)  for  the  rapid  pressure- strain  correlation  was  first  suggested  in  Girimaji  (2002)  by 
simply  truncating  (clipping)  the  LRR-QI  model.  Here,  by  testing  the  model  in  different  flows,  it 
was  confirmed  that  the  truncated  LRR  model  i  s  fully  realizable  for  all  initial  conditions  in  the 
Lumley  triangle.  In  the  near  future,  we  will  undertake  the  development  of  fully  realizable  and 
physically  superior  models  based  on  the  new  constraints. 
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Figure  1.6  Selected  trajectories  (es, 
model  at  S*  =  5. 


A-,  B,  and  s)  in  PS  flow  calculated  with  the  LRR-QI 
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Figure  1.8  Selected  trajectories  (es,  Bi,  CO,  a,  Jr,  and  »)  in  AC  flow  calculated  with  the  SSG 
model  at  S*  =  50. 


Figure  1.10  Contribution  of  three  parts  of  flie  Reynolds  stress  equation  in  trajectory  evolution 
shown  at  stages  of  Mya  and  RS  violations.  Case;  PS  flow  at  S*=5\  LRR-QI  model  describes  rapid 
pressure  strain  correlation. 


Figure  1.11  Contribution  of  three  parts  of  the  Rejoiolds  stress  equation  in  trajectory  evolution 
shown  at  stages  of  Mija  and  RS  violations.  Case:  AE  flow  at  S*=15;  LRR-QI  model  describes 
rapid  pressure  strain  correlation. 


Figure  1.16  Selected  trajectories  calculated  with  the  LRR-QI  model  in  AE  flow  at  S  =  50 
Denotations :  (*,  h)  initial  conditions,  (X)  equilibrium  point,  0)  locations  of  the 
A4w-constraint  violation,  (®)  RS-realizability  violation  point. 


Figure  1.18  Trajectories  obtained  with  the  original  LRR-QI  model  (— )  and  truncated  LRR 
model  ( — )  for  two  initial  conditions  (a)  and  (b)  in  PS  flow  at  o  =  5. 


Figure  1.19  Trajectories  obtained  with  die  original  LRR-QI  model  ( — )  and  truncated  LRR 
model  ( — )  for  two  initial  conditions  (a)  and  (b)  in  AE  flow  at  5  =15. 


Figure  1.20  Contribution  of  three  parts  of  the  Reynolds  stress  equation  in  evolution  of  original 
and  corrected  LRR  trajectories  shown  at  stages  of  Miju  and  RS  violations.  Case:  PS  flow  at 
S*=5. 


Figure  1.21  Contribution  of  three  parts  of  the  Reynolds  stress  equation  in  evolution  of  original 
and  corrected  LRR  trajectories  shown  at  stages  of  Muki  and  RS  violations.  Case:  AE  flow  at  S  = 


Chapter  2  Pressure-Strain  Correlation  in  Homogeneous  Anisotropic 
Turbulence  Subjected  to  Rapid  Strain-Dominated  Distortion 


2.1  INTRODUCTION 

Accurate  modeling  of  die  pressure- strain  correlation  continues  to  be  one  of  the  important 
challenges  in  the  second-moment  closure  mediod  of  turbulence  modeling.  In  particular,  the 
“rapid”  portion  of  the  pressure-strain  correlation  has  been  the  subject  of  many  analytical  and 
modeling  studies. 

When  turbulence  is  subjected  to  rapid  distortion  (RD),  the  description  of  turbulence  evolution 
simplifies  significantly  leading  to  the  so-called  rapid-distortion  theory  (RDT)  equations  which 
are  linear  in  fluctuating  velocity.  In  this  limit,  the  physics  of  rapid  pressure-strain  correlation  can 
be  studied  in  relative  isolation  as  the  complicating  effects  of  slow  pressure-strain  correlation  and 
dissipation  are  absent.  Another  important  feature  of  this  limit  is  that  the  fluid  behaves  as  an 
elastic  rather  than  viscous  medium.  The  stresses  and  other  turbulence  properties  including  the 
pressure-strain  correlation  depend  on  the  total  strain  (deformation)  experienced  by  the  fluid 
element  rather  than  the  current  strain  rate.  Beyond  this,  our  understanding  of  the  behavior  of 
initially  anisotropic  turbulence  subject  to  rapid  distortion  is  relatively  poor  as  much  of  the 
previous  RD  studies  have  focused  on  initially  isotropic  turbulence. 

The  focus  of  this  work  is  on  physics  of  pressure-strain  correlation  in  homogeneous  anisotropic 
turbulence  subject  to  rapid  strain-dominated  distortion.  Specifically,  we  study  extreme  cases  of 
initial  turbulence  anisotropy:  two-component  (2C)  and  one-component  (1C)  turbulence,  along 
with  isotropic  turbulence.  This  choice  is  motivated  by  the  argument  that  if  rapidly-distorted 
turbulence  behavior  can  be  well  understood  and  modeled  at  the  three  extremes  of  the  Lumley 
invariant  triangle  map’,  the  behavior  at  any  arbitrary  level  of  anisotropy  can  be  inferred  by 
suitable  interpolation.  Several  types  of  strain-dominated  rapid  deformation:  homogeneous  shear 
(HS),  plain  strain  (PS),  axisymmetric  contraction  (AC),  and  axisymmetric  expansion  (AE),  are 
considered  in  this  study.  Once,  some  degree  of  understanding  of  the  rapid  pressure-strain 
correlation  physics  is  developed,  we  will  attempt  to  extract  some  general  guidelines  to  aid  in 
future  model  developments. 

Rotation-dominated  (elliptic  flows)  mean  flows  and  flow  in  rotating  frames  are  examined  in  a 
separate  study  as  the  behavior  of  these  flows  in  flie  RD  limit  is  significantly  different  firom  that 
of  strain-dominated  flows. 

Modeling  Issues 

The  rapid  pressure- strain  correlation,  in  general,  depends  on  the  velocity  spectrum  tensor  and, 
hence,  cannot  beau  nique  f  unction  o  f  R  eynolds  s  tresses  a  nd  m  ean  v  elocity  g  radients* .  In  the 
RDT  calculations  (multi-point  closure)  pressure-strain  correlations  appears  in  closed  form  as  the 
velocity  spectrum  is  fully  known.  To  accurately  describe  pressure- strain  correlation  at  the  one- 
point  closure  level,  at  a  minimum  two  independent  turbulence  field  tensors  -  componentality  (C) 
and  dimensionality  (D)  —  are  needed’^’  Componentiality  refers  to  the  number  of  non- zero 


diagonal  components  in  die  Reynolds  stress  tensor,  whereas  dimensionality  is  a  statement  about 
the  wavenum^r  vector:  it  refers  to  the  number  of  orthogonal  directions  in  which  wave  numbers 
are  permitted'^.  Pressure-strain  correlation  can  possibly  take  a  large  range  of  values  for  a  given 
combination  of  Reynolds  stresses  and  mean  velocity  field,  depending  on  the  dimensionality 
tensor*^.  In  flie  traditional  second  moment  closure  (SMC)  approach,  this  all  important 
dimensionality  tensor  is  not  knoAvn.  Without  the  knowledge  of  dimensionality  tensor,  the 
pressure-strain  correlation  closure  problem  is  not  well-posed  as  the  solution  is  not  unique.  The 
challenge  of  traditional  one-point  closure  modeling  is  to  select  a  single  appropriate  value  within 
the  range  of  allowable  values  for  the  pressure- strain  correlation. 

Given  this  limitation,  no  single  one -point  closure  model  can  simulate  the  entire  range  of  physics 
incumbent  in  the  RDT  equations.  The  best  that  can  be  expected  of  any  traditional  one -point 
closure  model  is  that  it  captures  the  most  fundamental  features  of  the  RDT  equations.  As  a  result, 
in  one-point  closure  modeling,  many  of  die  pathological  aspects  of  the  RDT  behavior  will  have 
to  be  excluded  in  favor  of  die  more  dominant  features.  In  this  paper,  we  will  focus  only  on 
modeling  issues  pertaining  to  realizability  and  fundamental  consistency  with  rapid  distortion 
theory. 

Our  goal  is,  rather,  to  develop  modeling  constraints  that  go  beyond  the  Crow  constraint  and 
enjoy  a  wider  range  of  validity  and  applicability  in  rapidly  distorted  turbulence. 


Simple  analysis  shows  that  pressure- strain  correlation  can  be  written  as  following*  in  rapidly 
distorted  isotropic  turbulence: 


(2.1.1) 


where 
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Sif  =- 
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1  r  dUj  ^  dUj 

^ dxj  dxi 


Pij=-<  UiUk  >  — <  u  jUk  >  — • 

OXj^ 


In  the  above:  is  the  “rapid”  part  of  the  pressure-strain  correlation;  is  the  “rapid”  part 

of  the  pressure  fluctuation;  17, •  and  Ui  are  mean  and  fluctuating  velocities  components; 
k=^< UfUi  >  is  the  turbulent  kinetic  energy;  Py  is  the  production  tensor;  P  =  ^ Pa ;  and  dy  is 

the  Kronecker  delta  s5mibol.  Equation  (2.1.1)  is  called  the  Crow  constraint  and  its  use  in 
pressure-strain  correlation  modeling  can  be  traced  to  Rotta^. 


Currently,  coefficients  in  pressure-strain  correlation  models  are  chosen  to  satisfy  (2.1.1)  at  the 
RD  limit.  Models  satisfying  the  Crow  constraint  perform  well  in  rapidly-distorted  isotropic 
turbulence.  The  models,  however,  are  much  less  satisfactory  (unphysical  and  unrealizable)  in 
anisotropic  turbulence,  as  will  be  shown  later.  As  rapidly-distorted  anisotropic  turbulence  is  of 
great  practical  interest,  it  is  important  to  develop  modeling  constraints  such  as  equation  (2.1.1) 
for  anisotropic  turbulence. 


Previous  research 

A  review  of  previous  rapid-distortion  modeling  studies  is  given  in  Hunt  and  Carruthers^. 
Sreenivasan  and  Narasimha^  and  Maxey^  (and  references  therein)  investigated  the  influence  of 
various  types  of  distortion  and  anisotropy  levels  on  turbulence  evolution.  In  their  studies,  only 
special  cases  of  the  initial  energy  spectrum  tensor  were  considered.  Various  cases  of  three- 
component  (3C)  axisymmetric  turbulence,  with  the  Reynolds  stress  component  along  the  axis  of 
symmetry  increasing  from  its  isotropic  value  of  2/3A:  to  the  maximum  allowed  value  of  k,  were 
studied.  Maxej^  investigated,  in  some  detail,  the  role  of  the  pressure-strain  correlation  at  the 
rapid  distortion  limit  and  the  performance  of  the  model  of  Launder  et  al.^  in  die  homogeneous 
shear  flow.  Hunt  &  Carruthers^  further  pointed  out  that  the  asymptotic  value  of  Reynolds  stresses 
depend  on  the  initial  anisotropy  level.  An  important  conclusion  from  all  diese  studies  is  that  the 
level  of  initial  anisotropy  plays  a  crucial  role  in  the  subsequent  turbulence  evolution.  However, 
these  studies  stop  well  short  of  developing  general  modeling  constraints  or  guidelines  in  rapidly 
distorted  anisotropic  turbulence. 


Current  study 

The  objectives  of  the  present  study  are  to: 

1)  Develop  a  better  understanding  of  the  pressure- strain  correlation  in  highly  anisotropic 
turbulence  subject  to  various  types  of  rapid  distortion; 

2)  Assess  model  performance  (physical  consistency,  realizability)  by  comparison  with  RDT  data; 

3)  Seek  physical  princ^les  and  mathematical  constraints,  which  can  be  used  to  guide 
development  of  better  pressure- strain  correlation  models  in  the  RD  limit  for  anisotropic 
turbulence; 

4)  To  estimate  bounds  on  die  values  of  the  pressure-strain  correlation. 

Our  first  objective  is  similar  to  that  of  previous  studies  mentioned  above.  In  the  present  study,  we 
proceed  further  to  characterize  important  physical  features  of  pressure-strain  correlation  and 
establish  new  modeling  guidelines  in  anisotropic  turbulence.  These  guidelines  should  lead  to 
pressure- strain  correlation  models  of  high  degree  of  fidelity  to  the  Navier-Stokes  equations  in  the 
RD  limit. 


2.2  GOVERNING  EQUATIONS 

The  velocity  fluctuation  evolution  equations  in  the  RD  limit  are  die  following  (see,  e.g..  Pope*) 


Duj  _  dUj  1  d 
Dt  ’  dxi  p  dxj 


lv2pW 

p  dxj  dxi 


(2.2.1) 


Here,  The  equations  are  written  in  Cartesian  coordinates. 

Dt  '  Sxi  dt  ^ 

In  homogeneous  turbulence,  die  velocity  and  the  pressure  fields  can  be  written  in  terms 
of  their  Fourier  components: 


k 


(2.2.2) 


where  icQ)  is  the  wavenumber  vector  and  u(  ic.t )  is  the  Fourier  coefficient  vector  of  the  velocity 
fluctuation.  Then,  equations  (2.2.1)  transform  to  the  evolution  equations  for  the  components  of 

A 

vectors  K{t)  and  u{t) 


dKj  dU, 

— -  =  -^l — 
dt  dXj 


(2.2.3) 


dUi 

dxk 


subject  to  the  incompressibility  condition: 


(2.2.4) 


upCi^O.  (2.2.5) 

In  Cambon  and  Scott®,  it  is  suggested  that  equations  (2.2.3)-(2.2.4)  be  referred  to  as  the  Kelvim 
Townsend  equations.  For  given  initial  conditions,  equations  (2.2.23)-(2.2.5)  can  be  solved 
directly.  Then,  the  covariance  of  two  Fourier  coefficients 

kij(Jc,t)  =<u*  (K,t)ujiK,t)  >  (2.2.6) 


can  be  extracted  from  the  data  for  each  given  K(t) .  Summation  of  (2.2.6)  over  all  wavenumber 
vectors  gives  the  Reynolds  stress  components  in  the  physical  space: 


Anotiier  way  to  obtain  ^,y  and  Reynolds  stresses  is  to  solve  directly  flie  evolution  equation  for 
the  covariance  of  two  Fourier  coefficients  of  a  given  wavenumber  vector 


dRjj 

dt 


dx^  ^  dx^  dxk  k-2 


+  ^Rjk 


dUi  KjKi 


(2.2.8) 


which  is  derived  using  equation  (2.2.4).  Details  of  the  derivation  can  be  found  in  Pope®  (p.  412), 
for  instance.  Now,  the  continuity  equation  takes  the  following  form 


KiRij(K)=KjRij(K)  =  0. 


(2.2.9) 


Equation  ^.2.8)  is  a  numerically  more  efficient  alternative  to  equation  (2.2.4)  for  calculating 
Reynolds  stresses  at  the  RD  limit*®  and,  hence,  has  been  adopted  in  this  study.  Summation  of 
(2.2.8)  over  all  icQ)  yields  the  Reynolds  stress  evolution  equation  in  the  RD  limit: 


d  < UfUj  > 
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=  P  +  0) 
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where  production  Py  is  given  by 


( .  dUj  .  dUi  ^ 


(2.2.10) 


and  die  “rapid”  pressure-strain  correlation  ^  is 


=  2S 


Rik 


K  V 


dUi  KjKi 


+  Rjk 


dUj  tCjKi^ 


Thus,  the  RDT  equations  can  be  used  directly  to  study  the  physics  of  production  and  pressure- 
strain  redistribution  processes. 

In  the  current  study,  RDT  data  are  used  to  evaluate  directly  the  accuracy  of  three  popular  models 
for  the  “rapid”  part  of  the  pressure-strain  correlation  -  IP  (isotropizatiomof-production)  model**, 
LRR  modef,  and  SSG  modeP  -  for  given  isotropic  and  anisotropic  initial  conditions.  All  models 
can  be  represented  by  the  general  form: 


(2.2.11) 


=  -C\Pb,j^{fl-C\^)kS,j+C,k\b,,Sj,  -jbMj 

+  C,k^^W„*b„W„) 

In  the  above,  by  denotes  the  anisotropy  tensor  and  Wy  is  flie  mean  rotatioitrate  (vorticity) 


tensor: 
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dUj  dU  j 
dxj 


The  models  differ  only  in  the  value  of  coefficients,  which  are  given  in  Table  2.1.  The  value  0.8 
of  the  coefficient  C?  recommended  in  all  three  models,  comes  from  constraint  (2.1.1). 


2.3  CALCULATION  PROCEDURE 

The  RDT  equations  ((2.2.3),  (2.2.8)-(2.2,9))  and  modeled  Reynolds  stress  evolution  equation 
(2.2.10)  are  solved  in  a  variety  of  homogeneous  turbulent  flows.  Mean  flows  investigated  here 
are: 


homogeneous  shear  (HS) 


dxj 


"0  S  0^ 
0  0  0 
^0  0  0; 


plain  strain  (PS) 


0  0^ 
0-5  0 
^0  0  0; 


axisymmetric  contraction  (AC) 


dU, 

dXj 


0  0  ^ 

0  -s/2  0  , 

^0  0  -5/2; 


and  axisymmetric  expansion  (AE), 


dU, 

dXj 


^-25  0  0^ 
0  5  0 

^  0  0  5^ 


Initial  condition  specification  requires  the  definition  of  terms  componentiality  (C)  and 
dimensionality  (D)  of  the  turbulence  field.  In  the  present  paper,  turbulence  dimensionality  is 
always  equal  to  three  unless  prohibited  by  incompressibility  condition. 

In  keeping  with  the  main  objective  of  the  paper,  highly  anisotropic  as  well  as  isotropic  turbulent 
velocity  fields  are  chosen  as  initial  conditions.  Anisotropic  cases  considered  are:  one-component 
(1C)  and  axisymmetric  two-component  (2C)  turbulent  states.  Three  kinds  of  1C  initial  conditions 
are  considered: 


They  are  denoted  as  ICl,  1C2,  and  1C3  respectively.  The  three  axisymmetric  2C  initial 
conditions  are 


which  are  denoted  as  2C1, 2C2,  and  2C3  respectively.  Isotropic  condition  is  given  by 

by- =0,/,  7  =1,2,3.  (2.3.3) 

Values  of  lie  anisotropy  tensor  given  in  (2.3.1)-(2.3.3)  are  used  directly  to  solve  modeled 
Reynolds  stress  evolution  equation  (2.2.10). 

To  solve  RDT  equations  0.2.3),  (2.2.8),  and  (2.2.9),  initial  turbulence  fields  with  desirable 
properties  should  be  generated.  Since  many  different  choices  of  the  wavenumber  vectors  and  the 
Fourier  coefficients  of  the  velocity  fluctuation  can  yield  flie  required  anisotropy,  the  choice  of 
the  initial  fields  for  these  quantities  must  be  made  carefully.  In  the  present  study  the  following 

assumptions  on  tiie  initial  distributions  of  vectors  K(t  =  (i)  and  u(K,t  =0^are  made: 

a)  Velocity  fluctuations  in  all  permissible  directions  are  equally  energetic.  That  is,  Fourier 
coefficients  of  initial  velocity  fluctuations  in  all  permissible  directions  have  the  equal 
magnitude.  Whether  Fourier  coefficient  vectors  are  permitted  in  a  specific  direction  is 
determined  fi'om  the  initial  turbulence  componentiality. 

b)  For  a  given  vector  u(t  =  0)  all  permissible  wavenumber  vector  directions  are  equally 
probable.  Permissible  wavenumber  vector  directions  are  determined  purely  by 
incompressibility  condition  (2.2.5). 


Thus,  the  most  general  (unbiased)  initial  velocity  field  can  be  generated  for  the  specified 
Reynolds  stress  tensor.  For  1C  turbulence,  the  initial  vector  u  is  aligned  along  the  corresponding 

A 

unit  vector.  In  2C  turbulence,  initial  vectors  u  are  uniformly  distributed  in  a  circle  on  the 
permissible  plane.  In  3C  i  sotropic  case,  vectors  u  are  uniformly  distributed  on  a  sphere.  The 
magnitudes  of  u  are  such  that  the  total  initial  turbulent  kinetic  energy  is  recovered. 

For  each  vector  u(t  =  0)  ,iihe  corresponding  wavenumber  vectors  K(t  -0)  are  determined  from 
assumption  (b)  above.  Because  equation  ^.2.8)  does  not  depend  on  the  magnitude  of  the 
wavenumber  vector,  but  on  its  direction  only,  one  can  assume,  without  loss  of  the  generality  that 
all  permissible  wavenumber  vectors  are  of  equal  magnitude  initially.  Thus,  to  generate  initially 
isotropic  turbulence,  wavenumber  vectors  are  distributed  evenly  on  the  surface  of  a  unit  sphere. 
In  1C  turbulence,  vectors  if  (t)  are  uniformly  distributed  in  the  unit  circle  lying  in  the  plane 

normal  to  the  initial  vector  m  .  In  2C  turbulence,  the  generation  of  vectors  ic  is  explained  below 
for  the  case  of  2C3  turbulence. 

In  2C3  turbulence  initial  vectors  «  are  equally  distributed  in  the  plane  (1,2)  (Fig.  1).  The 
corresponding  wavenumber  vector  distribution  is  dictated  by  incompressibility  condition 
(assumption  b).  Effectively,  the  full  set  of  vectors  ic  is  divided  into  families,  with  each  family 

having  the  same  value  of  the  projection  of  ic  on  axis  3  (Fig.  1):  kTj  =|r:|cos(0  '~^+dd),  for 

i=2,N  and  JC3=|ff|,for  i  =  1.  Here,|fr|  =  l,  0*=O,and  =;r/2(A^-l^.  All  families  cany 
equal  amount  of  energy.  Then,  knowing  the  number  of  vectors  ic  in  each  family,  one  can 
determine  the  energy  associated  with  each  vector  ic  and  Rjj(ic,t  =  0).  Vectors  ic  in  the  cases  of 
2C2  turbulence  and  2C1  turbulence  are  generated  in  the  similar  manner. 

A  fourth-order  Runge-Kutta  scheme  was  used  for  time  integration  of  all  equations. 


2.4  RESULTS  AND  DISCUSSION 

Results  firom  computations  of  RDT  equations  ((2.2.3),  (2.2.8)-(2.1.9))  and  modeled  Reynolds 
stress  evolution  equation  (2.2.10)  are  now  presented.  The  mean  flows,  the  initial  conditions,  and 
basic  outcomes  from  the  RDT  and  model  calculations  are  given  in  Table  2.2.  The  RDT 
calculations  show  evolution  of  anisotropy  level  in  some  cases  (denoted  by  ‘E’)  and  no  evolution 
in  others  (‘N’).  All  isotropic  turbulence  cases  evolve  from  their  initial  states. 


Our  RDT  data  interrogation  focuses  on  the  role  of  the  pressure-strain  correlation,  especially  its 
relation  to  production.  The  analysis  is  centered  around  the  anisotropy  evolution  equation  at  the 
RD  limit: 
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dt 
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where  P'\j  (=Pij-  - 
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8yP)  is  the  production  of  anisotropy.  Also  a  crucial  value  in  the 


analysis  is  the  magnitude  of  anisotropy,  which  evolves  according  to 


1  ^byhy  1  ^  ('>1)  ) 

2  dt  2k^ 


(2.4.2) 


We  first  investigate  if  a  generalized  form  of  constraint  (1)  is  possible  for  an  anisotropic  flow,  that 
is,  if  the  pressure-strain  correlation  merely  counteracts  certain  fraction  of  the  production.  In  such 
a  case,  the  pressure-strain  correlation  could  have  one  of  two  forms: 


where  the  proportionality  coefficient  in  each  expression  could  be  a  constant  or  a  scalar  invariant 
of  the  mean  velocity  gradient  tensor.  The  simplest  maimer  to  verify  the  existence  of  such  a  linear 
relationship  would  be  to  compute  the  following  ratios 


(r) 


R\  = 


P"  P" 
^  ij^  IJ 


or  R2  = 


<D 

ptpt 


(2.4.3) 


The  first  ratio  /?1  in  (2.4.3)  is,  probably,  more  justified  in  anisotropic  turbulence  considering  tiie 
form  of  anisotropy  evolution  equation  (2.4.1)  in  the  RD  limit.  In  isotropic  turbulence,  we  have 
/?l  =  i?2  =  0.36  in  accordance  with  (1).  The  RDT  data  confirm  that  flie  ratio  i?l  is  indeed 
preferable  to  R2  because  at  t  =  0  Rl  s  hows  1  ess  v  ariation  f or  d  ifferent  i  nitial  flow  c  onditions. 
Nevertheless,  even  the  ratio  Rl  changes  firom  case  to  case:i?l  =  l  for  HS  (1C2,  2C1)  and 
Rl  =  0.25  for  HS  (2C2),  PS  (2C1,  2C2,  2C3),  AE  and  AC  (2C2,  2C3)  flows.  Further  Rl  and  m 
change  substantially  with  turbulence  evolution.  Some  examples  of  Rl  evolution  are  shown  in 
Fig.  2.2.  In  none  of  the  cases  considered,  the  latio  Rl  (and  R2)  is  equal  to  the  isotropic  value 
0.36. 


One  clear  conclusion  can  be  drawn  from  RDT  data:  the  pressure- strain  correlation  does  more 
than  merely  counteract  the  anisotropy  production.  In  tiie  absence  of  a  unique  relationship 
between  the  production  and  pressure-strain  correlation,  subtler  connections  are  now  investigated. 
One  of  the  main  features  of  the  pressure-strain  correlation  is  its  redistributive  nature.  The  fact 

that  the  pressure-strain  correlation  termO,^'^'^'’  is  traceless  reveals  that  it  does  not  alter  the  total 

energy.  It  is  generally  believed  tiiat  tiie  pressure-strain  correlation  removes  energy  from  high- 
energy  components  and  enhances  lower-energy  components.  In  contrast,  the  production  process 
is  believed  to  increase  anisotropy  by  injecting  energy  into  selected  components  of  the  Reynolds 
stress  tensor.  The  specific  issues  investigated  are  the  following. 

1)  Can  there  be  redistribution  without  production?  That  is,  what  happens  to  when 

R'',  =  0. 


2)  Does  the  “rapid”  pressure- strain  correlation  always  extract  energy  from  high-energy 
components  and  deposit  it  into  low-energy  ones?  This  can  be  answered  by  monitoring  the 

sign  of  (refer  to  equation  2.4.2). 

3)  Can  more  energy  be  redistributed  than  produced  in  the  first  place?  That  is,  can  the  /?1  (or 
Rl)  value  be  greater  than  unity. 

2.4.1  RDT  results 

Non-evolving  cases 

The  RDT  calculations  demonstrate  that  in  some  cases  (denoted  by  ‘N’  in  Table  2.2)  turbulence 
anisotropy  does  not  change  from  its  initial  level.  In  other  words,  the  term  dby  / dt  on  the  left 
side  of  equation  (2.4.1)  is  equal  to  zero  for  each  component  of  the  anisotropy  tensor  b,j.  It  is 
straightforward  to  show  (see  the  Appendix)  that  in  such  cases  all  components  of  the  production 
term  P"ij  are  exactly  equal  to  zero.  It  follows  ttien  from  equation  (2.4.1),  ftiat  tiie  pressure-strain 

tensor  components  should  also  be  equal  to  zero.  Thus,  if  P",y  =  0,  then  =  0.  In  order  to 
confirm  this,  we  examine  the  second  invariant  of  tensors  and  ,  ie.,  P",jP''y  and 

.  The  second  invariant  contains  information  about  all  tensor  components.  If  all 
individual  components  of  the  tensor  are  equal  to  zero,  so  is  the  invariant.  The  RDT  calculations 
show  that  when  the  invariant  P”,jP”y  is  equal  to  zero,  flienO/'^O),/'^'’  is  also  equal  to  zero. 
Note  that  in  these  cases,  I)j  or  P',y  need  not  vanish. 


Evolving  cases 

Cases  in  which  ftie  RDT  predicts  turbulence  anisotropy  evolution  are  denoted  by  ‘  E’  in  T able 
2.2.  If  the  pressure-strain  correlation  decreases  anisotropy  and  the  production  term  increases  its 

level,  flien  in  (2.4.2)  should  be  negative,  and  P''ijby  should  be  positive.  The  computed 

RDT  values  of  P"^fblJ  and  ^J'^by  in  some  representative  cases  are  plotted  in  Fig.  3.  The 

P  'ijby  evolution  is  shown  with  a  solid  line  (without  symbols),  while  the  ^J'^by  evolution  is 

represented  by  dashed  line  (wifliout  symbols).  It  is  clearly  seen  that  P"yby  is  always  positive  as 

expected,  while  <i>J'^by  is  predominantly  negative  (again  as  expected).  In  one  case  (Fig.  3e)  it 

assumes  small  positive  values  at  latter  stages  of  simulation.  It  can  be  argued  though  that  at  this 
late  stage  of  evolution,  RDT  equations  may  not  represent  the  real  turbulence  physics.  It  would  be 

interesting  to  analyze  the  circumstances  under  which  ^^'^by  is  positive.  As  a  general  rule,  we 
conclude  that  <^y'^by  is  negative,  implying  that  the  pressure-strain  correlation  predominantly 
reduces  flie  Reynolds  stress  anisotropy. 

The  RDT  data  demonstrate  that  the  magnitude  of  Rl  andR2  initially  never  exceeds  unity.  Some 
results  are  shown  in  Fig.  2.  In  few  cases,  for  example  in  AE  (2C2)  flow  (Fig.  2b),  Rl  exceeds 
unity  at  long  times.  This  is  again  at  time  when  the  RDT  may  not  accurately  reflect  the  turbulence 


physics.  The  conclusion  then  is  that  the  magnitude  of  the  redistribution  does  not  in  general 
exceed  tiiat  of  the  production. 

For  strain-dominated  mean  flows,  RDT  results  appear  to  imply  that  the  behavior  of  the  pressure- 
strain  correlation  is  as  follows  irrespective  of  the  initial  anisotropy: 

1)  4)/"'' =0, when  P'V  =  0; 

2) ®/%  <  0; 

3)  .Rl  andi?2  <1. 

It  would  be  reasonable  to  demand  that  models  be,  at  least,  qualitatively  consistent  with  die  RDT 
data. 

2.4.2  Model  calculations 

Three  models  for  the  “rapid”  part  of  the  pressure-strain  correlation  (P,  LRR,  and  SSG)  are  now 
compared  with  he  RDT  data.  The  most  fundamental  aspects  of  the  model  calculations  are 
tabulated  in  Table  2.2.  In  the  table,  ‘ER’  stands  for  an  evolving  realizable  trajectory  and  ‘EU’ 
denotes  an  evolving  unrealizable  solution.  The  table  demonstrates  that  in  many  cases  models 
yield  evolution  of  turbulence  anisotropy  when  the  RDT  clearly  shows  no  evolution.  In  other 
cases,  when  the  RDT  shows  evolution,  model  behavior  is  umealizable.  Thus,  two  categories  of 
undesirable  model  behavior  can  be  recognized:  inconsistency  with  turbulence  (RDT)  physics  and 
unrealizable  evolution  trajectories. 

Inconsistency.  If  a  turbulence  model  predicts  evolution  of  Reynolds  stress  anisotropy  when  RDT 
indicates  no  evolution,  the  model  is  labeled  as  being  inconsistent  with  turbulence  physics.  In 
Table  2.2,  when  a  model  shows  evolution  (denoted  as  ‘ER’  and  ‘EU’),  but  the  RDT  does  not 
(‘N’),  we  have  inconsistency. 

Unrealizability.  As  is  common  practice,  a  model  is  labeled  as  unrealizable,  if  it  produces 
negative  diagonal  components  of  the  Re)molds  stress  tensor  (in  the  principal  coordinates).  In 
such  a  case,  a  model  evolution  trajectory  passes  from  the  realizable  region  (inside  Lumley 
triangle)  to  the  non-realizable  area  (outside  Lumley  triangle).  In  all  cases  identified  with  ‘EU’  (in 
Table  2.2),  the  model  behavior  is  unrealizable. 

We  will  now  compare  model  calculations  with  RDT  data  and  try  to  explain  the  physical  failing 
underlying  inconsistent  and  unrealizable  model  behavior.  In  particular,  the  connection  between 
flie  violations  and  the  proposed  modeling  guidelines  will  be  sought. 


Non-evolving  cases 

We  discuss  in  detail,  at  first,  the  results  from  four  cases:  AC(1C2),  AE(1C2),  HS(lCl),  and 
PS(lCl).  The  cases  are  identified  by  the  mean  flow  and  the  initial  turbulence  anisotropy  (in 
parentheses).  In  fliese  cases,  RDT  calculations  show  no  evolution,  as  is  seen  from  Table  2.2. 

Thus,  bij(t)=bij(0);  F'ijbij=0  and  <Pi/’’^bij  =  0 .  Models,  however,  demonstrate  various 
behavior  patterns  including  inconsistent  and  unrealizable  evolution.  The  model  evolution 


trajectories  are  tracked  on  the  (|,  77)  invariant  plane:  where  the  invariants  are  defined  by 

677  =  byb ji ,  6|  =  bybjf^  bjfi.  To  assess  realizability  violation,  the  Lumley  invariant  triangle  is 

also  shown.  A  trajectory  is  deemed  unrealizable  when  it  crosses  the  boundary  of  the  Lumley 
triangle. 

Case  AC  (1C2).  All  models  predict  spurious  evolution  in  this  case,  but  evolution  trajectories  stay 
inside  the  Lumley  triangle  (Fig.  2.4a).  Thus,  all  three  models  considered  predict  inconsistent,  but 
realizable  behavior.  To  understand  the  origin  of  poor  model  behavior,  we  now  examine  the  ratio 

i?l  (Fig.  2.4b)  and  the  invariants  /*”,y6,yand  <I),y^''^6,y(Fig.  2.4c).  The  initial  value  of 
implied  by  each  model  is  non- zero,  while  the  F'yF'y  value  is  equal  to  zero  as 
can  be  seen  from  i?l  00  at  t  =0  (Fig.  2.4b).  Non-zero  model  frien,  causes  the  model 

trajectory  to  evolve.  In  RDT  calculations,  both  and  F'y  are  equal  to  zero  resulting  in  no 

evolution.  The  cormection  between  inconsistent  model  behavior  and  violation  of  criterion  1  is 
easily  e  stablished.  F  rom  F  ig.  2.4(c),  it  is  seen,  that  for  each  model,  F'yby  is  positive  and 

is  negative.  To  summarize,  all  models  violate  criteria  1  and  3,  and  not  criterion  2  in 

this  case. 

Case  HS  (ICl).  The  IP  model  predicts  no  anisotropy  evolution  in  agreement  with  the  RDT.  In 
both  cases  -  IP- model  and  RDT  calculations  -  initial  values  of  and  F'y  are  equal  to  zero. 

Results  from  LRR  and  SSG  models  are  shown  in  Fig.  2.5.  The  LRR  model  shows  inconsisteit 
and  unrealizable  behavior.  The  SSG  model  predictions  are  inconsistent,  but  realizable  (Fig. 
2.5a).  From  Fig.  5(b)  it  is  clear  that  both  models  violate  criteria  1  and  3.  As  regards  criterion  2, 
flie  LRR  model  violates  it,  but  the  SSG  model  does  not  (Fig.  5c).  The  inference  is  that 
unrealizable  model  behavior  is  linked  to  violation  of  criterion  2.  This  conclusion  is  further 
corroborated  by  following  cases. 

Case  PS  (ICl).  In  this  case  all  models  exhibit  inconsistent  and  unrealizable  behavior  (Fig.  2.6a). 
All  three  models  violate  all  three  criteria  (Figs.  2.6b  and  2.6c). 

CaseAE  (1C2).  Again,  all  models  predict  inconsistent  and  unrealizable  behavior  (Fig.  2.7a),  and 
all  three  criteria  are  violated  (Figs.  2.7b  and  2.7c).  The  SSG  model  further  predicts  that 
production  tends  to  reduce  anisotropy  {F’y  by  <0),  which  is  completely  inconsistent  with  the 

RDT  and  turbulence  physics. 

Evolving  cases 

Now  we  consider  cases  where  the  RDT  shows  evolution  of  initial  turbulence  anisotropy  (see 
Table  2.2).  In  all  these  cases,  all  three  models  show  anisotropy  evolution  also.  In  this  sense  tiieir 
behavior  is  consistent  with  the  RDT.  When  a  model  predicts  consistent  and  realizable  behavior, 
no  criterion  violation  is  observed.  Some  examples  are  shown  in  Figs.  2.3  (a-c)  for  criterion  2  and 
in  Figs.  2.2(a)  and  2.2(b)  for  criterion  3.  Results  for  all  cases  studied  are  summarized  in  Table 


2.3.  In  many  cases,  the  model  evolution  is  unrealizable.  As  examples,  cases  AC  (2C2)  and  PS 
(2C2)  will  be  discussed  in  detail.  (Criterion  1  does  not  apply  for  evolving  cases). 

Case  AC  (2C2).  The  RDT  trajectory  starts  from  the  2C-axisymmetric  turbulence  comer  of  the 
Lumley  triangle  and  evolves  towards  the  1C  turbulence  comer  (Fig.  2.8a).  All  three  model 
trajectories  initially  exit  the  Lumley  triangle,  but  later  recover  and  re-enter  the  realizability  zone 
(Fig.  2.8a).  The  model  prediction  after  the  first  trajectory  exit  from  the  Lumley  triangle  is 
irrelevant  and,  therefore,  the  model  behavior  must  be  deemed  unrealizable.  The  reason  for  the 

violation  can  be  seen  in  Fig.  2.3(d).  The  value  of  predicted  by  each  model  is  initially 

positive  leading  to  an  increase  of  anisotropy  level.  At  later  times,  the  values  become 

negative  as  the  Reynolds  stresses  recover  to  realizable  values.  The  RDT  value  of  <bij^^by  is 
always  negative.  Criterion  2  is  violated,  but  not  criterion  3  (Fig.  2.2c). 

Case  PS  (2C2).  The  RDT  trajectory  is  similar  to  the  AC(2C2)  case  (Fig.  2.8b).  All  models  are 
again  unrealizable  (Fig.  2.8b).  The  degree  of  violation  in  this  case  is  much  larger  and  there  is  no 
recovery  to  realizable  values.  The  explanation  is  again  found  by  investigating  the  behavior  of 

<bij'^^by  presented  in  (Fig.  2.3e).  The  model  values  of  ^ij'^^by  are  positive  and  substantially 

larger  than  in  the  previous  case  making  recovery  to  realizable  values  difficult.  The  RDT  values 
are,  again,  always  negative.  Criterion  3  is  not  violated  (Fig.  2.2d). 

All  the  calculated  results  are  summarized  in  Table  2.3.  In  Table  2.3,  ‘Cl’,  ‘C2’,  and  ‘C3’ 
represent  the  three  criteria.  Violation  of  a  criterion  is  denoted  by  ‘V’  and  norv  violation  is  shown 
by  ‘N’.  To  clearly  demonstrate  the  coruiection  between  undesirable  behavior  and  criteria 
violation,  two  columns  -  recording  consistency  and  realizability  -  are  also  added  in  Table  2.3.  In 
these  columns,  again  ‘V’  indicates  violation  and  ‘N’  shows  norvviolation.  The  relationship 
between  undesirable  model  behavior  and  modeling  guideline  violation  is  clearfy  evident. 
Violations  of  criteria  1  and  3  are  associated  with  spurious  evolution  and  inconsistent  model 
behavior.  Violations  of  criterion  2  and  realizability  condition  are  connected. 

Implications  of  the  RDT  and  model  calculations  are  now  summarized. 

1)  Without  exceptions,  flie  RDT  data  indicates  that  =0,  when  P''y  =  0.  Inconsistent 

model  behavior  is  imminent  when  is  non-zero  as  F’y  vanishes.  Therefore,  criterion  1 

is  a  mathematical  constraint  that  a  model  must  be  mandated  to  satisfy.  Another  important 
inference  is  that  for  this  combination  of  Reynolds  stresses  and  mean  velocity  gradients,  the 
“rapid”  pressure-strain  correlation  assumes  a  unique  value,  which  is  zero. 

2)  From  flie  results  in  this  section,  it  is  clear  that  Oy^^^by  being  positive  is  not  a  desirable 
feature  in  a  model,  since  all  observed  realizability  violations  are  associated  with  such 
behavior.  In  all  these  cases,  flie  RDT  values  of  Oy^'^^by  remain  negative.  The  underlying 

physics  is  now  explained.  When  <by^''^by  is  positive,  the  model  extracts  energy  from  low- 
energy  components  and  deposits  it  into  high-energy  components.  This  action  continues  even 


when  the  low-energy  component  is  completely  depleted,  driving  that  component  into 
negative  values.  Mathematically,  this  is  much  like  counter- gradient  transport,  which  almost 
always  has  a  destabilizing  influence.  While  the  RDT  data  does  permit  this  destabilizing 
energy  transfer  in  very  rare  instances,  physics  incumbent  in  the  RDT  equations  ensures  a 
delicate  balance  between  production  and  redistributive  processes  leading  to  realizable 
behavior.  It  remains  to  be  seen  if  single-point  closure  models  can  be  developed  to  replicate 
this  delicate  balance.  Until  this  issue  is  addressed  satisfactorily,  criterion  2  must  be 
considered  as  an  important  constraint  that  a  model  should  satisfy. 

3)  Much  like  criterion  2,  criterion  3  appears  to  be  physically  reasonable,  although  not  always 
trae  in  RDT.  This  criterion  also  suggests  approximate  bounds  on  the  pressure-strain 

correlation  magnitude:  0  <  Thus,  the  larger  the  anisotropy 

production,  the  wider  is  the  range  of  values  diat  the  correlation  can  take.  The  margin  of  error 
incurred  in  pressure-strain  correlation  modeling  is  therefore  likely  to  increase  with  increasing 
anisotropy  production,  even  with  the  most  accurate  models. 

2.5  CONCLUSIONS 

Rapid  distortion  calculations  of  initially  anisotropic  turbulence  are  performed  to  better 
imderstand  the  physics  of  the  pressure-strain  correlation.  The  specific  objectives  are  to  establish 
“universal”  principles  that  underlie  pressure-strain  correlation  behavior  and  estimate  the  extent 
of  range  of  values  that  the  correlation  can  take  for  given  combinations  of  Reynolds  stresses  and 
mean  velocity  gradients. 

In  the  rapid  distortion  limit,  turbulence  evolution  depends  on  various  factors:  type  of  anisotropy, 
degree  of  anisotropy,  and  the  imposed  mean- flow  deformation.  Consequently,  the  pressure-strain 
correlation  term  will  also  depend  on  these  parameters.  While  the  behavior  of  the  pressure-strain 
correlation  can  be  quite  complex  and  not  easily  generalized  in  anisotropic  turbulence,  certain 
physical  features  can  be  identified  as  being  universal.  The  RDT  calculations  performed  in  this 
paper  show  that  in  mean  strain-dominated  flows  the  dynamics  of  turbulence  is  such  that: 

i)  <I)/"^=0,when  P'V  =  0. 

ii) P"yby>  0,  whereas  <  0. 

iii)  0  ^  ^  F’ijF'ij. 

It  is  proposed  fliat  the  above  criteria  be  used  as  guidelines  for  developing  “rapid”  pressure-strain 
correlation  models. 

Reynolds  stress  transport  model  calculations  performed  at  the  rapid  distortion  limit  display 
unphysical  (inconsistent)  and  unrealizable  behavior  for  some  initial  turbulence  anisotropies  and 
mean  flows.  Inconsistent  behavior  refers  to  spurious  evolution  predicted  by  a  model  when  the 
RDT  indicates  no  change  in  anisotropy  from  its  initial  level.  Clearly,  in  a  full-turbulence 
equation  set,  there  may  be  evolution  from  an  anisotropic  initial  state.  But  this  evolution  must  be 
initiated  by  “slow”  pressure-strain  correlation  or  dissipation  terms. 


Close  examination  of  the  RSTM  results  reveal  ttiat  both  unrealizable  and  inconsistent  model 
behavior  can  be  directly  related  to  violations  of  the  constraints  proposed  in  this  paper. 
Inconsistent  behavior  occurs  when  the  model  pressure-strain  correlation  does  not  vanish  as  the 
production  of  anisotropy  goes  to  zero.  Unrealizable  model  behavior  is  observed  when  the  model 
pressure-strain  correlation  unphysically  attempts  to  increase  anisotropy  rather  dian  diminish  it. 

In  ongoing  work''*  we  demonstrate  that  current  models  yield  realizable  and  consistent  behavior 
when  their  coefficients  are  modified  to  satisfy  the  constraints  proposed  in  this  paper. 

The  results  and  analysis  presented  in  this  paper  should  aid  in  better  understanding  die  role  of  the 
pressure-strain  correlation  in  turbulence  and  hence,  enable  us  to  develop  improved  models  for 
this  term.  Similar  examination  of  mean  rotation-dominated  flows  will  be  undertaken  in  the 
future. 
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APPENDIX: 


The  RDT  results  show  (Table  2.2)  that  in  the  following  cases:  ICl,  1C3  (for  all  mean  flows), 
1C2  (AC,  AE,  PS),  2C1(AC,  AE),  and  2C2(HS),  initial  anisotropy  of  turbulence  does  not 
change.  These  are  cases,  when  the  production-ofanisotropy  term  P"y  in  equation  (2.4.1)  is 

equal  to  zero  for  all  Rejnolds  stresses  at  t  =  0 . 


As  an  example,  let  us  consider  the  mean  flow  under  axisymmetric  contraction,  where  i/j  ,  =  5, 
^1,1  =^3,3  =-5/2,  and  the  rest  of  the  components  is  equal  to  zero.  In  such  a  flow,  the 
components  of  the  Reynolds  stress  production  tensor 

„  du, 

Pij=-<  >  - - <  UjUt  >  — 

dx^  dx^ 

are  Pj ,  =  -2S  <«,«!>,  P22  =  5  <  M2M2  > »  P33  =  5  <  M3M3  > ,  and  the  others  are  equal  to  zero. 
For  different  initial  turbulence  anisotropy  cases  one  has 
ICl:  < MjMi  >=2k,  P,i  =-4Sk,  P22  =^33  =  0,  2P  =  -4Sk, 

1C2:  <  U2U2  >=  2k ,  P22  =  2Sk ,  P, j  =  P33  =  0,  2P  =  2Sk, 

1C3:  <  1/31/3  >=2k,  P33  =  2Sk,  /ji  =^22  2P  =  2Sk, 

2C1:  <t/,Mj  >=0,  <  U2U2  >=<M3M3  >=k,  Pji  =0,  P22  =P33  =  Sk,  2P  =  2Sk, 

2C2:  <  U2U2  >=  0 ,  <  M,w,  >=< «3«3  >=k.  Pi  1  =  —2kS,  P22  =0,  P2^=Sk ,  2P  =  -Sk , 

2C3:  <  >=  0 ,  <  i/ii/j  >=<  U2U2  >=k.  Pi  1  =  -2kS,  P22  =  ,  P33  =  0 ,  2P  =  -Sk . 

2 

Then,  the  components  of  the  production-of-anisotropy  tensor  P"^-  =  P^-  -2Pbij-—5yP  in  each 
case  are 

ICl:  P"„  =  P”22  =  P"33  =  0, 

1C2:  P"„  =  P"22  =  P"33  =  0, 

1C3:  P", ,  =  P''2,  =  P''33  =  0, 

2C1:  P”„  =  P"22  =  P"33=0. 

2C2:  F\^  =  -^Sk,  P"22  =  0,  P"i^=^Sk, 

2C3:  P”,i  =  -|5A:,  P’'22  =  ^Sk,  P"33  =  0. 

Non-diagonal  components  of  P'^  are  equal  to  zero  in  all  cases.  It  is  seen  that  only  in  two  cases 
(2C2  and  2C3)  there  are  non- zero  components  of  P",y  and,  therefore,  the  invariant  P”,y  P",y  has 
a  non- zero  value.  In  these  two  cases,  the  RDT  data  show  the  turbulence  anisotropy  evolution 


from  the  initial  level.  In  other  cases  (ICl,  1C2,  1C3,  and  2C1)  there  is  no  production  of 
anisotropy  and,  as  RDT  computations  demonstrate,  no  change  of  initial  anisotropy  occurs.  1 
follows  then  from  equation  (2.4.1),  that  the  pressure-strain  tensor  components  have  to  be  also 
equal  to  zero  in  such  cases. 

For  other  mean  flow  configurations,  the  similar  analysis  can  be  easily  made,  which  leads  to  the 
same  conclusions. 
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Table  2.1  Model  coefficients  for  expression  (12) 
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Table  2.2  Notations:  (E)  evolving  trajectory;  (N)  no  evolving;  (U)  unrealizable  evolution;  (R) 
realizable  evolution.  1C3  (AC,  AE)  case  is  equivalent  to  1C2  (AC,  AE);  and  2C3  (AC,  AE)  case 
is  equivalent  to  2C2  (AC,  AE). 
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Table  2.3  Notations:  (Cl)  criterion  1;  (C2)  criterion  2;  (C2)  criterion  3  (C3); 


(N)  non- violation;  (V)  violation  . 


Figure  2.1  Coordinate  scheme  for  the  2C3  initial  turbulence  state. 

Figure  2.2  Evolution  of  the  ratio  Rl.  Cases:  a)  HS  (1C2),  b)  AE  (2C2),  c)  AC  (2C2),  d)  PS 
(2C2). 

Notations  are  (— )  RDT,  (#')  IP  model,  (SD)  LRR  model,  («)SSG  model. 

Figure  2.3  Evolution  of  terms  F'ij  by  and  ^i/^^by  in  equation  (17).  Cases:  a)  PS  (isotropic), 
b)  HS  (1C2),  c)  AE  (2C2),  d)  AC  (2C2),  e)  PS  (2C2).  Notations:  for  F'y  by  (— )  RDT,  (O — ) 

IP,  (A—)  LRR,  (D— )  SSG;  for  (-  -)  RDT,  (0-)  IP,  (A-)  LRR,  (□-)  SSG. 

Figure  2.4  AC  flow  with  1C2  initial  turbulence  state:  a)  Lumley  triangle,  b)  evolution  of  the  Rl 
parameter,  c)  evolution  of  invariants  P"y  by  and  ^J’^^by .  In  Fig.  4(a)  the  black  dots  (•)  mark 

the  comers  of  Lumley  triangle  corresponding  to  isotropic,  2C,  and  1C  conditions.  Other 
notations  see  in  Figs.  2.2  and  2.3. 

Figure  2.5  HS  flow  widi  ICl  initial  turbulence  state:  a)  Lumley  triangle,  b)  evolution  of  the  Rl 
parameter,  c)  evolution  of  invariants  F'y  by  and  (notations  are  the  same  as  in  Figs. 

2.2-2.4). 

Figure  2.6  PS  flow  with  ICl  initial  turbulence  state:  a)  Lumley  triangle,  b)  evolution  of  the  Rl 
parameter,  c)  evolution  of  invariants  F'y  by  and  (notations  are  the  same  as  in  Figs. 

2.2-2.4). 

Figure  2.7  AE  flow  with  1C2  initial  turbulence  state:  a)  Lumley  triangle,  b)  evolution  of  the  Rl 
parameter,  c)  evolution  of  invariants  F'y  by  and  (notations  are  the  same  as  in  Figs. 

2.2-2.4). 

Figure  2.8  Solution  trajectories  in  the  Lumley  triangle  Cases:  a)  AC  (2C2  and  isotropic),  b)  PS 
(2C2  and  isotropic).  Notations  are  (— )  RDT,  (O)  IP  model,  (A)  LRR  model,  (□)  SSG  model  for 
initially  anisotropic  turbulence.  The  comers  of  Lumley  triangle  are  marked  with  (•). 


Fig.  2.1 


3 


1 


Fig.  2.3 


CHAPTER  4 


4.1.  Turbulent  Transport  Models  for  PANS 


The  two-equation  PANS  model  as  derived  by  Girimaji,  2003  [1]  is  as  follows 


„  dK, 
dt 
deu 
dt 


dxi 


(1) 


+  (2) 

where  the  two  terms  for  transport  of  unresolved  kinetic  energy  and  dissipation,  Tk^  and  are 


Tfcu 


dK^ 


dxj  Wfc  dxj )  ^  ^  Qxj 

dxjXa^dxj)  ^  ^  dxj 


(3) 

(4) 


In  the  above  equations,  the  unclosed  terms  that  require  modeling  are  the  convection  of  Ku,  f-u  by 
resolved  fluctuations  {Uj  —  Uj).  This  is  considered  next.  The  equations  for  7*^,  and  can  be 
written  wholly  in  terms  of  PANS  variables  by  expressing  the  ratio  of  total  viscosity  to  turbulent 
Prandtl  numbers  as 

t/t  ^  ^  1  f,  ^  vuf.  ,5. 

(<^fc,o-e)  (crfc,a«)  e„  fl  {crk,(Tt)fk 

In  order  to  close  these  transport  terms  completely,  we  examine  two  extremes  or  limits. 
Zero-Transport  assumption:  This  assumption  is  based  on  the  fact  that  the  resolved  fluctuations 
{Uj  —  Uj)  do  not  take  part  in  the  net  transport  of  unresolved  kinetic  energy  (Ku)  at  the  scales  of 
resolution  in  which  we  are  interested  (see  Fig  1.). 

Let,  Vj  =  {Uj  —  Uj),  then 


VjKu  = 

VjCu  — 


A  and 
B 


(6) 

(7) 


—  dx 

'5(-T> 

„  dx 
Ku{x-—) 

^u{x-  —) 


““  dx 
Vjix+-) 

dx 

Kuix+-) 

^uiX+^) 


Figure  1:  Control  Volume  showing  flux  of  and  €« 


Now,  the  first  assumption  implies  {A,  B)  «  constant  in  the  scales  of  interest  and  hence, 


d{A,B) 


=  0  in  the  large  scales  of  interest 


Thus,  we  are  only  left  with 


which  is  finally  written  as 


T  =  ^  ff'ufcdKu\ 

dxj  Vcrjfe/I  dxj  ) 

T  =  ^ 

dxjKa^fldxj) 


T  = 

dxj  Vo-fc„  dxj  j 

T  = 

dXj  Va£„  dxj 


with  the  unresolved  Prandtl  numbers  expressed  as 


=  7-o-fc 
/« 

=  4‘^e  (14) 

Je 

Msiximum-Transport  Model:  The  second  case  is  based  on  a  more  generalized  gradient-transport 
closure  and  states  that  the  resolved  field  velocity  is  proportional  to  the  resolved  field  eddy  viscosity 
Ur,  which  can  be  defined  as 

Ur  =  Ut—Uu  =  Cy, - Cy - 


(15) 


Thus, 


^  dKu  __  {yt  -  Vu)  dKu 

(Tk  dxj  Ck  dxj 

(16) 

Ur  deu  iUt-Uu)d€u 

(Tf  dXj  <7e  dXj 

(17) 

and  hence  ,  ,  >  „ 

d{A)  _  d  r{ut-  Vu)  dKu\ 

dxj  dxj  V  (Tk  dxj  J 

(18) 

d{B)  _  d  ^{vt-Vu)d€u\ 
dxj  dxj  \  cTt  dxjj 

(19) 

From  equations  (3.3)  and  (3.4),  we  get 

.  a  /  vt  d{Ku,eu)\  d  r{vt-Vu)d(Ku,eu)\ 

\  tu)  dxj\{(Tk,(Te.)  dxj  )  dxj\  {ak,(Tt)  dxj  J 

_  d  /  Vu  9{Ku)€u)\ 

dxj\{ak,(Te)  dXj  J 

(20) 

which  is  expressed  explicitly  as. 

^  d  f  Vu  dKu\  J 

=  dxjWdxj) 

(21) 

T  = 

dxj\at^dxjj 

(22) 

with  the  modified  Prandtl  numbers  for  unresolved  kinetic  energy  and  dissipation  as, 


CTku  =  <^k  (23) 

(24) 

If  we  consider  a  continuum  (Fig.l)  whose  width  (da;)  corresponds  to  an  order  of  length  equal  to 
the  intermediate  scales  of  motion,  then  the  flux  of  (or  eu/cu)  between  regions  1  and  2  of  the 
continuum  is  given  by 

4"  )Ku{x  +  — )  —  Vj(x  —  )Ku(x  — )  (25) 

Expanding  each  term  in  a  Taylor  series  and  after  cancelling  various  terms,  we  are  left  with  the 
convection  term  for  K^,  and  similarly  1^-^.  Now  since  e„  is  a  small-scale  quantity  and 

consequently  its  statistics  converge  to  a  mean  value  quickly,  for  a  given  width  of  the  continuum,  it  is 
reasonable  to  assume  that  its  mean  changes  slowly  between  regions  1  and  2,  and  hence  its  gradient 
(or  transport)  is  neghgibly  small  in  the  region  of  interest.i.e. 

dtu 
dxj 


0 


(26) 


(27) 


Thus,  in  equation  (3.4)  the  only  remaining  term  is  the  viscous  transport  term 

r 

dxj\(Ttdxj) 

where  the  ratio  of  the  total  viscosity  (i/j)  to  the  dissipation  Prandtl  number  (cr^)  can  be  written  as 


,,  r  ^  C  ^  V 

a  =  bLj-  =  (28) 

which  is  essentially  the  Zero-TVansport  model  for  unresolved  dissipation 

In  the  equation  for  Ku  however,  transport  terms  due  to  both  viscosity  and  resolved  fluctuations 
exist,  as  Ku  is  a  large  scale  quantity  whose  mean  changes  considerably  between  regions  1  and  2  and 
therefore  its  gradient  (or  transport)  between  regions  1  and  2  of  the  continuum  is  significant  and 
cannot  be  neglected.  Hence,  in  this  case  we  have 


Tk 


dxj 


(29) 


and  therefore  make  use  of  the  second  assumption  based  on  gradient-transport  closure  and  obtain 


g  /  i/f  dKu\ _ d  ^ {ut  -  dKu\ 

dxj  VcTfc  dxj  )  dxj  \  cTfc  dxj  J 

d  /i/udKu\ 
dxj  VcTfc  dxj  / 


(30) 


which  is  essentially  the  Mziximum-Transport  model  for  unresolved  kinetic  energy,  Ku-  This 
proves  that  the  ZT  model  is  better  suited  for  modeling  unresolved  dissipation  (c„),  while  the  MT 
model  is  adequate  for  unresolved  kinetic  energy  (Ku). 

Objective:  The  objective  of  the  current  study  is  to  determine  which  of  these  two  transport  models 
performs  better  for  Ku  and 


4.2.  RESULTS  AND  DISCUSSION 


4.2.1.  Maximum-TVansport  Vs  Zero-TVansport  Model 

In  this  section,  the  imsteady  and  turbulent  flow  results  from  numerical  simulations  conducted  on 
the  three-dimensional  lid-driven  cavity  flow  for  a  spanwise  aspect  ratio  (SAR  =  2L/W)  of  3:1:1  are 
presented  and  compared  with  published  experimental  observations  and  LES  (Jordan,  1994  [2]).  The 
results  present  comparisons  made  between  simulations  conducted  using  the  Maximum-TVansport 
(MT)  and  the  Zero-TVansport  (ZT)  models  to  ascertain  which  of  the  models  performs  better  for  Ku 
and  c„.  Qualitative  comparisons  are  available  in  terms  of  contours  of  X-vorticity,  while  quantitative 
comparisons  with  LES  are  realized  through  plots  of  Centerline  Mean  Velocities  at  a  z  =  0.28  plane, 
although  experimental  results  are  available  at  the  symmetry  plane. 


All  figures  have  been  generated  on  completion  of  the  flow  simulation  corresponding  to  a  flow  time 
(T)  of  approximately  115  seconds.  The  plots  for  statistical  quantities  along  cavity  centerlines  have 
been  generated  such  that  for  the  first  10  cycles  (corresponding  to  about  60  seconds  of  flow  time), 
no  sampling  of  data  is  done  in  order  for  the  flow  to  become  fully  turbulent,  after  which  data  is 
sampled  every  6.5  seconds  to  obtain  realizations  of  mean  quantities.  Results  have  been  presented 
for  fk  values  of  1.0  (RANS),  0.7  and  0.4. 

Figures  2  and  3  show  snapshots  of  the  mean  U  and  V  velocity  profiles  along  the  cavity  centerlines, 
while  Fig.  4  presents  contours  of  X-vorticity  for  a  typical  RANS  (/fc=l)  simulation.  Very  few  flow 
features  are  visible  in  the  contour  plots  which  illustrates  the  fact  that  RANS  is  inherently  incapable 
of  capturing  many  of  the  flow  features  and  the  results  for  statistical  quantities  show  its  inability  to 
produce  accurate  results,  especially  near  the  boundaries  of  the  cavity. 


Mean  U-Velocity 


Figure  2:  Mean  U-Velocity  profiles,  bottom  wall  closeup,  (RANS,  UG) 


Mean  V-Velocity 


Center  Line  Mean  Velocity  Profiles — RANS, Uniform  Grid 


Figure  3:  Mean  V-Velocity  profiles,  downstream  wall  closeup,  (RANS,  UG) 


Figure  4:  X-Vorticity  Contours  along  downstream  wall,  (RANS,  UG) 


The  next  set  of  figures  (5-13)  show  comparisons  of  the  two  turbulent  transport  models,  ZT  and 
MT,  with  each  other  as  weU  as  against  LES  and  experiment.  Prom  the  statistical  comparisons  for 
=  0.7  (Figs.  5,  6),  its  appears  that  (i)  both  models  seem  to  be  performing  better  than  RANS 
and  show  greater  convergence  of  the  PANS  calculations  towards  LES  results;  and  (ii)  the  plots  also 
illustrate  the  fact  that  the  ZT  model  performs  marginally  better  than  the  corresponding  MT  model 
near  the  cavity  boundaries.  Although  marginally  better  accuracy  is  observed  from  the  statistical 
comparisons,  instantaneous  X-Vorticity  contours  (Figs.  7  and  8)  along  the  downstream  wall  display 
the  fact  that  considerably  more  scales  of  flow  are  resolved  by  the  ZT  model  when  compared  to  the 
MT  model,  while  both  models  are  able  to  uncover  significeintly  more  scales  of  flow  than  RANS,  in 
terms  of  appearance  of  Taylor  Goertler  (TGL)  vortex  pairs.  About  three  pairs  of  TGL  vortices  are 
observed  for  the  ZT  model  (Fig.  7)  while  only  one  pair  can  be  seen  for  the  MT  model  (Fig.  8),  for 
a  value  of  fk=0.7. 
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Center  Line  Mean  Velocity  Profiles-ZTM, Uniform  Grid 


Figure  5:  Mean  U- Velocity  profiles,  bottom  wall  closeup,  fk  =  0.7,  (ZT  Vs  MT) 
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Center  Line  Mean  Velocity  Profiles— ZTM  Vs  MTM, Uniform  grid 


Figure  6:  Mean  V-Velocity  profiles,  downstream  wall  closeup,  /*,  =  0.7,  (ZT  Vs  MT) 


As  the  value  of  fk  is  reduced,  more  and  more  pairs  of  vortices  are  uncovered  and  the  number  of  flow 
features  resolved  by  the  ZT  model  is  considerably  more  than  the  corresponding  MT  model  for  a 
given  fk,  speciflcally,  5  pairs  of  TGL  vortices  are  observed  for  the  ZT  model  as  against  3  pairs  for  the 
MT  model,  for  fk=0A  (Figs.  11,  12).  Statistical  comparisons  for  the  models  indicate  that  although 
both  models  seem  to  perform  fairly  well  along  the  center  of  the  cavity,  the  ZT  model  appears  to 
perform  better  than  the  MT  model  in  the  layers  adjoining  the  upper  and  lower  boundaries.  Further, 
as  the  value  of  fk  is  reduced,  the  ZT  model  shows  closer  agreement  with  LES  and  experiment  than 
the  corresponding  MT  model,  for  a  given  fk-  This  trend  is  evident  in  flgure  9,  10  {fk  =  0.4). 


Center  Line  Mean  Velocity  Profiles — CTM, Uniform  Grid 


Figure  9:  Mean  U- Velocity  profiles,  bottom  wall  closeup,  fk  =  0.4,  (ZT  Vs  MT) 
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Center  Line  Mean  Velocity  Profiles — 2TM  Vs  MTM, Uniform  grid 


Figure  10:  Mean  V-Velocity  profiles,  downstream  wall  closeup,  fk  =  0.4,  (ZT  Vs  MT) 


Figure  12:  X-Vorticity  Contours  along  downstream  wall,  /*  =  0.4,  (MT,  UG) 


Finally,  two  important  inferences  that  can  be  drawn  from  theses  results  are  (i)  the  ZT  model 
performs  better  than  the  corresponding  MT  model  for  a  given  fk  and  (ii)  as  the  value  of  fk  is 
decreased,  better  convergence  towards  LES  and  experiment  is  observed.  The  second  inference  is 
especially  evident  in  plots  13-16.  Figures  17-20  show  a  comparison  between  plots  of  contours  of 
X-vorticity  for  various  fk  values  and  essentially  display  the  fact  that  a  reduction  in  the  value  of  fk 
causes  more  and  more  scales  to  be  resolved. 


Center  Line  Mean  Velocity  Profiles — ^ZTM, Uniform  Grid 


Figure  13:  Mean  U- Velocity  profiles,  bottom  wall  closeup,  (ZT,  UG) 
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Figure  15;  X-Vorticity  Contours:  fk  =  0.2 


Figure  16:  X-Vorticity  Contours:  fk  =  0.4 


Figure  17:  X-Vorticity  Contours:  fk  =  0.7  Figure  18:  X-Vorticity  Contours:  RANS 

4.3.  CONCLUSIONS  AND  SUMMARY 


In  the  present  study,  a  detailed  numerical  investigation  of  the  three-dimensional  lid-driven  cavity 
flow  is  carried  out  using  a  new  turbulence  modeling  technique,  know  as  the  Partially  Averaged 
Navier  Stokes  method  or  PANS.  Flow  studies  have  been  conducted  at  Re=10000  on  a  lid-driven 
cavity  for  a  spanwise  aspect  ratio  (SAR)  of  3:1:1,  with  the  top  surface  imparted  with  a  constant 
velocity  of  1  m/s.  The  investigation  primarily  focuses  on  the  effects  of  two  turbulent  transport 
models  (the  Zero-Transport  (ZT)  and  Maximum-IVansport  (MT))  on  PANS  calculations  to  ascer¬ 
tain  which  model  performs  better  for  and  e„. 


The  main  conclusions  that  can  be  drawn  by  observing  all  the  results  are:  (i)  as  the  value  of  fk  is 
decreased  from  1  (RANS)  to  0  (DNS)  more  scales  of  flow  are  resolved  in  terms  of  the  appearance  of 
flow  features,  such  as  the  TGL  vortices,  as  evident  from  the  contour  plots  of  X-vorticity,  and  more 
importantly  (ii)  the  ZT  model  appears  to  perform  better  than  the  MT  model.  As  a  final  note,  in 
view  of  all  the  results,  it  is  seems  justifiable  to  conclude  that  the  PANS  model  is  capable  of  serving 
as  abridging  model  between  RANS  and  DNS. 
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